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1 Introduction 

Most of the available information about matter in 
our Universe, and in our Galaxy in particular, comes 
indirectly from the collection of the electromagnetic 
radiation (from meter waves to 7 rays) that was emit- 
ted or absorbed by this matter. A completely differ- 
ent information is provided by the cosmic ray nuclei, 
which constitute a genuine sample of galactic matter. 
Many different nuclei species are observed, in a wide 
range of energy, and with different origins. Some 
of them come unaltered from the sources (they are 
called primaries), otliers (secondaries) come from 
nuclear reactions between the primaries and the in- 
terstellar medium, or from the disintegration of un- 
stable species. Moreover, the trajectories of these 
nuclei from creation to detection are rather erratic, 
due to the influence of the galactic magnetic field on 



sible to follow tlie direction of an incoming nucleus 
back to the source. If we were able to understand 
clearly the processes by which all these nuclei are 
produced, accelerated and propagated in the Galaxy, 
the wealth of data available now or in the near fu- 
ture would yield most valuable information about the 
matter content and magneto-hydrodynamical prop- 
erties of our Galaxy. In principle, it would even be 
possible to discover some evidence for new physics 
(e.g. supersymmetry) or new objects (e.g. primor- 
dial black lioles or stars made of antimatter) as they 
can give rise to the emission of charged antinuclei 
and make an extra contribution to the observed cos- 
mic ray fluxes. This review presents a summary of 
the work made in this direction by the authors from 
1999 to 2002, in a lapth-isn-iap collaboration. As 
a first step, we tried to reach a quantitative under- 
standing of the propagation of cosmic ray nuclei in 
the energy range 100 MeV/nuc-100 GeV/nuc. More 
precisely, we described propagation with a diffusion 
model, in which the free parameters are adjusted to 
account for the available data on cosmic rays. This 
provides the regions of the parameter space allowed 
by the data. As a second step, we took advantage 
of this model to investigate several points concerning 
astrophysics and astroparticle physics. 

During this study, various aspects related to the 
"standard" or to more speculative processes were ex- 
amined in detail. These may be summarized as fol- 
lows: 

• Standard cosmic rays 

— Diffusion parameters from secondary-to- 
primary ratio 

— The flux of standard secondary antiprotons 
and antideuterons 

^ Spatial origin 

— Radioactive species and the local bubble 

— Evolution of composition with energy 

• Exotic cosmic rays 

— Baryonic Dark Matter 

— Antimatter 

— Supersymmetric particles 

— Primordial Black Holes 

The extraction of the diffusion parameters is the cen- 
tral goal since all conclusions follow from their values. 



2 Propagation models 
2.1 The context 

The study of cosmic ray radiation raises a great num- 
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charged nuclei are accelerated and how do they prop- 
agate in chaotic magnetic fields. These questions 
are actually still actively studied and debated, but 
the nature of the sources and propagation media on 
which cosmic ray physicist focus has slightly changed 
in the last twenty years. The pioneering work of 
Parker Q, who first used the diffusion-convection 
equation, focused on Solar modulation. The same 
equation was used in the late seventies to in- 
vestigate the acceleration of Galactic Cosmic Rays 
(hereafter GCR, in opposition to Solar Cosmic Rays, 
SCR) in sources that were strongly suspected to be 
supernovae (SN) remnants. Whereas SCR deal with 
the Sun magnetic field and plasma in the Solar cav- 
ity, GCR have to face the question of transport in 
the largely unknown galactic magnetic field. These 
studies reached an even larger scale, since the dis- 
covery of ultra high energy cosmic rays (UHECR) 
that are certainly extragalactic in origin. The ex- 
act nature of this radiation is not clearly established 
but if these are charged nuclei, the relevant medium 
in which these nuclei propagate is the extra-galactic 
magnetic field. From the UHECR source point of 
view, one has to imagine some powerful astrophysical 
sources where "standard" acceleration occurs, but an 
astroparticle solution (i.e. heavy meta-stable decay- 
ing new particles) could furnish the energetic parti- 
cles as well. New particles could thus be discovered 
by the study of cosmic radiation, which is an inter- 
esting coming back to the very first concern of this 
field of research, when the muon and the positron 
were first observed in cosmic rays. 

Notwithstanding the fact that UHECR is one of 
the driving subjects in the development of the so- 
called astroparticle field, the study of much less en- 
ergetic particles, with energies ranging from GeV to 
PeV, may also put constraints on the existence of 
new particles. As discussed at length in the follow- 
ing, these low energy particles have a galactic ori- 
gin, and by their study, many questions concern- 
ing the Galaxy may be addressed, such as LiBeB 
primordial abundances, dark matter content of our 
Galaxy, or the nature of the sources. It might be 
possible to reach a consistent picture of either con- 
ventional (SNs, Wolf-Rayet stars) or possible more 
exotic (micro-quasars anti-globular clusters [Q) 
galactic sources. As a consequence, we underline 
that GCR nuclei in the GeV-PeV energy range pro- 
vide a quite interesting laboratory for both astro- 
physicists and particle physicists. From now on, we 
will only consider cosmic rays in this energy range. 

The great amount of data in various energy ranges 
has led to a quite good understanding of charged nu- 
clei propagation (see e.g. j^) along with the induced 
7-ray production at low and high latitude in the 
galactic plane. As the gyration radius of charged par- 
ticles in the galactic magnetic field is small, the prop- 
agation is intimately related to the detailed structure 
of this turbulent magnetic field. The latter is not 
observed directlv and one would like to use the cos- 



mic rays to infer its properties. It turns out to be a 
difficult task as, despite the success outlined above, 
several unknowns and inconsistencies remain at the 
quantitative level. Unlike the Solar case for which we 
have in situ observations of this turbulence - as early 
as in the mid-sixties (see e.g. |^) - favoring a Kol- 
mogorov spectrum, some recent MHD simulations 
along with our secondary/primary studies point to- 
wards greater spectral index of turbulence. Even 
though a satisfying global picture emerges, some nu- 
clei resist to a simple interpretation. If the first en- 
during problem ~ the depletion of the grammage dis- 
tribution at low values for sub-Fe/Fe ||, ^ ^ seems 
now to be solved thanks to new cross section mea- 
surements [Q, another lasting discussion is related 
to acceleration and selection mechanisms that lead to 
the observed abundances. Is it a chemical selection, 
i.e. First Ionization Potential bias (see e.g. ||ll|]), 
or a volatility bias related to grain destruction in 
SN explosions (12)? The accuracy of current data is 
not sufficient to conclude and the question remains 
in suspense, albeit the importance to elucidate the 
injection mechanism in acceleration models. It is 
at least known to a great certainty from radioac- 
tive primaries (^^Co, ^^Ni) synthesized in SN that 
10^ years have past between synthesis and accel- 
eration [|3|. The energy spectrum produced by the 
sources has also some indeterminacy. The accelera- 
tion models agree about a power-law dependence in 
rigidity TZ~°' (where the rigidity is defined as mo- 
mentum per unit charge TZ = p/Z) but numerical 
estimates of the spectral index a can be rather dif- 
ferent. Values a < 2.0 are preferred by accelera- 
tion theory (see and in particular ||l^ for a short 
and readable introduction on diffusive shock accel- 
eration), but if Kolmogorov spectrum for diffusion 
is retained, one is left with a « 2.5 from observed 
spectra that seems more problematic. Spectra may 
also differ from pure power laws. This depends also 
on the kind of sources involved (SN, explosion in 
wind bubbles, superbubbles) along with their plasma 
and magnetic states. There could be less explosive 
sources such a Wolf-Rayet stars which eject species 
through a powerful wind (this is the ^^Ne abundance 
anomaly, see e.g. [|l6|). 

This list is far from being exhaustive, even if we 
restrict ourselves to the case of ^GeV/nuc charged 
nuclei. To conclude on what we call "standard" CR 
studies, we emphasize that the limiting factor is re- 
lated to our poor knowledge of most production cross 
sections. 



2.2 Physical motivations for diffusive propa- 
gation 

The history of cosmic ray propagation theory traces 
back to the pioneering work of Fermi, which provides 
a statistical description of the way random mag- 
netic irregularities scatter the particles and acceler- 



2.3 Overview of the effects affecting propagation 
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the same decade, Chandrasekliar had shown rig- 
orously that diffusion could be equivalently described 
in terms of random walk, i.e. sequence of small er- 
ratic steps. By then, it was possible to imagine that 
the scattering mentioned above could lead to spa- 
tial diffusion. Actually, this hypothesis was first pro- 
posed on a phenomenological basis, as emphasized 



in Berezinskii et al. |f8 , in the context of Solar par- 
ticles transport, ft was then applied to cosmic ray 
transport in the Galaxy. The notion of diffusive mo- 
tion was confirmed and refined by more fundamental 
approaches, based on relativistic Boltzmann equa- 
tions [Q. ft now appears as a most valid description, 
that accounts for a wide range of observations (both 
in acceleration and propagation mechanisms). The 
reader is referred to |^ for a historical background 
about transport equations (see also |^ for a deriva- 
tion from the kinetic equation level and for the 
diffusion/convection equation). 

The linearized kinetic theory approach provides 
grounds for a consistent derivation of the transport 
equation, which reads, neglecting spallations and en- 
ergy losses for the sake of clarity 



|-V(A-V/- 



- VJ) - 



V.K 1 d 



p2 dp dp dp 



3 p^ dp 
dp 



(1) 



where / = f{t,r,p) is the phase space distribution. 
This equation contains most of the effects described 
below, like spatial diffusion, galactic wind, with the 
adiabatic energy loss associated, and diffusion in mo- 
mentum space. In this equation, spatial diffusion has 
been assumed to be isotropic. Actually, the diffusion 
coefficient K should be replaced by a tensor, with 
parallel and transverse components. As regards the 
first one, there is a strong consensus about a form 



i^ll(7^) = /vo/37^ 



2-K 



(2) 



where k is the spectral index of the turbulence spec- 
trum. The transverse component is still debated, and 
the two main propositions (given by quasi-linear or 
the Bohm conjecture) are probably wrong |Q. All 
the results presented here are based on the usual as- 
sumption that diffusion is isotropic, with a diffusion 
coefficient 



K{n) ^Kapn^ 



(3) 



where the normalization A'o and the spectral index 5 
should ideally be related to the astrophysical prop- 
erties of the interstellar medium. Unfortunately, our 
knowledge in this field is still demanding, and the 
value of the two parameters Kq and 5 can only be 
determined indirectly by the analysis of cosmic ray 
observations. 

The different diffusion schemes mentioned above 



of the reacceleration term Kpp The reader is 

referred to [|^, and we will not discuss this point 
further. 



2.3 Overview of the effects affecting propaga- 
tion 

This section is devoted to a brief overview of the 
physical effects that play a role in the propagation 
of cosmic rays. The next section (see will enter 
into more details and focus on the modelling of these 
effects. 



2.3.1 Geometry and content of the Galaxy 

The propagation of cosmic rays ceases to be of dif- 
fusive nature beyond some surface where they can 
freely stream out of the diffusive volume. The den- 
sity then drops to nearly zero, so that this surface 
may be considered as an absorbing boundary. The 
exact shape and dimensions of this boundary arc not 
known, but direct observations of the radio halo of 
external galaxies suggest that it might radially fol- 
low the galactic disc, with a greater thickness. In this 
work, the diffusive halo will be modelled as a cylinder 
of radius R = 2Q kpc and half-height L whose numer- 
ical value, still to be determined, is probably greater 
than a few kpc. The boundary thus imposes that the 
density satisfies N{r = R, z) = N{r, z = ±L) = 0. 
Embedded in this diffusive halo lies the galactic disc 
(hereafter "the disc" , in short) containing the stars 
and the gas . The gas is mostly made of hydrogen 
(90%), neutral and ionized, and hehum (10%) (the 
heavier nuclei that may be present are of negligible 
importance for our concerns). The different compo- 
nents (stars and gas) have different half heights hi, 
of the order of /i ~ 100 pc; they all satisfy h L, 
so that the disc will be considered as infinitely thin 
for all practical purposes. The density of interstellar 
matter is observed to be about tt-ism ~ 1 part cm~^ 
for all radii, so that we take n{r,z) = 2hd{z)nisM- 
Sources and interactions with matter are confined to 
the thin disc and diffusion which occurs throughout 
disc and halo with the same strength, is independent 
of space coordinates. The Solar System is located 
in the galactic disc (z = 0) and at a Galactocentric 
distance Rq ~ 8 kpc [^. A schematic view of the 
galactic model is shown in Fig. ^. 

2.3.2 Spallations: the importance of cross sections 

When a cosmic ray crosses the disc, it may interact 
with an interstellar hydrogen or helium nucleus and 
initiate a nuclear reaction (spallation). The impor- 
tance of this effect is governed by the corresponding 
cross sections. Actually, it is important to know not 
only the reaction or total cross section, which deter- 
mines the rate of destruction of a given CR species, 
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nel, which gives the formation rate of new nuclei. 
The cross section for a given channel is often referred 
to as the spallation or fragmentation cross section. 
The determination of all these cross sections is a nu- 
clear physics problem that shall be addressed in de- 



tails in Sec. 2.4.5 



2.3.3 Energy losses from interaction with the ISM 

There are two types of energy losses which are rele- 
vant for nuclei: ionization losses in the ISM neutral 
matter (90% H and 10% He), and Coulomb energy 
losses in a completely ionized plasma, dominated by 
scattering off the thermal electrons. The other ef- 
fects like brcmsstrahlung, synchrotron radiation and 
inverse Compton are negligible in the conditions con- 
sidered here. 

The Coulomb energy loss rate is given in pO, M 



dE 
'dt 



-iTTr^cnieC^ Z'^Ue In A 



Coul 



where 



In A w - In 
2 



In these expressions re and rue denote the classical 
radius and rest mass of the electron (Particle Data 
Group []), (ne) - 0.033 cm'^ and - lO'' K de- 
note the density and temperature of the interstellar 
electrons Z and M are the charge and mass 

numbers of the incoming nucleus and In A ~ 40 — 50 
is the Coulomb logarithm. 

The rclativistic expression giving the ionization 
losses is 



/ Ion 



27rr^mpc3z2 



s=H.Hc 



where 



S, EEE In 



^max — 



2TOeC^/3^7^Q 

72 



1 + [2 7 me /A/] 



2/32 



and /3oc ~ 0.01 c is the typical velocity of bound elec- 
trons in the hydrogen atom, Ig is the geometrical 
mean of all the ionization and excitation potentials of 
the considered atom, (/h = 19 eV and /hc = 44 eV), 
M ^ me is the incident nucleon mass, and is the 
density of the target atom in the ISM. 



2.3.4 Adiabatic losses from convective wind 

Among the phenomena affecting the propagation of 
cosmic rays, the magnitude of those presented in this 
section and the following are more subject to debate 
because of the uncertainty associated to parameters 
Vc and Vo- 
lt is very likely that the medium responsible for 
diffusion is moving away from the disc, with a veloc- 
ity Vc- This is referred to as convective or galactic 
wind, in analogy with the Solar wind. One of the 
effects of this galactic wind is to dilute the energy 
of the particle located in the disc in a larger vol- 
ume [ p9[ . This adiabatic expansion results in a third 
type of energy loss, depending on V.V^. Throughout 
our works, a very simple and tractable form for Vc 
is adopted, following |Q. It is considered to be per- 
pendicular to the disc plane and to have a constant 
magnitude throughout the diffusive volume, so that 
dVc/dz = except at z — where a discontinuity 
occurs, due to the opposite sign of the wind velocity 
above and below the galactic plane. In this case, the 
dilution effect on the energy is given by a term that 
can be expressed in the same form as ionization and 
coulomb losses 



J Adiab 



2m + ^fc \ V 
m + Ek J 3h 



(4) 



Ek stands for the total kinetic energy and it should 
not be confused with the kinetic energy per nucleus 
frequently used in this paper. We emphasize that 
this term corresponds to a process occurring only in 
the disc but not in the halo. 



2.3.5 Reacceleration 

Along with a the spatial diffusion coefficient K, 
Eq. (|^) also contains the momentum diffusion coeffi- 
cient Kpp] both are related to the diffusive nature of 
the process. The latter coefficient Kpp is related to 
the velocity of disturbances in the hydrodynamical 
plasma, called Alfven velocity. From the quasi-linear 
theory, it is given by (see e.g. 



'3(5(4- J2)(4_5) 



VaV/K{E).{5) 



In this expression, /iroac stands for the half-height of 
the cylinder in which reacceleration occurs. In our 
model, /ircac = hi but as Kpp depends only on the 
combination V^h-cca.c/h, the same diffusion in mo- 
mentum space is obtained for /ircac 7^ h, provided 
that the true Alfven velocity is given as a function 
of the parameter Va by 14 x h/hrcac- 



2.3.6 Full propagation equation 

The transport equation (p]) can be rewritten for each 
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N^{E) = dn^ /dE. As the momentum distribution 
function is normalized to the total cosmic ray num- 
ber density {n = 4tt J dp p^f), we have N^{E) = 
{in/ f3)p'^ to finally obtain, assuming steady-state 
(see Sec. |2.4.6D, 



-V 



KVN^E) - VcN^{E) 



dE 



^N^{E) 



_d_ 

dE 



-Otot 



= Q'iE) + 
dNi{E) 



dE 



where the following notation has been used for the 
total energy loss term 6tot = &ioss + &roac, with 



biossiE) 



dE 
'dt 



dE 



dE\ 



Ion \ / Coul V / Adiab 

and the reacceleration drift term defined as 



6rcac(i?) = ^^^jf ^ K, 



PP ■ 



We also use a compact notation to describe the most 
general form for a source term 

mfc>mj 

Q'{E)=q^Q^{E)+ r'^^iv^o), (7) 

k 




E, (GeV/rjc) 



Figure 1 : Characteristic times of several processes af- 
fecting the propagation of cosmic rays are displayed 
in the 100 MeV/nuc-100 GeV/nuc energy range. 
Typical values = 0.03 kpc^ Myr-\ b = 0.6 and 
Vc = 10 km s~^ were considered. The dominant pro- 
cess at energies higher than a few GeV is the escape 
through the boundaries of the diffusive volume. The 
effect of spallations is seen to be small for the prop- 
agation of protons, whereas it is crucial for heavy 
nuclei such as Fe. 



which includes primary sources - normalized abun- 
dance 50, spectrum Q^{E) -, but also s econd ary 
sources, coming from spallations (see Sec. 2.4.6 ) or 
radioactive decay of a heavier species (see Sec. 2.4.8 ). 
The relative magnitude of all the effects affect- 
ing propagation can be estimated from the typical 
timescale associated with these effects, as displayed 
in Fig. 0. 

Taking advantage of cylindrical symmetry and 
adding radioactive contributions localized in the disc 
and the halo, the previous equation may be rewritten 
as (making implicit the energy dependence) 



+2hS{z) {qlQ{E)q{r) - V {E)N' {r,Q) 



mk>mj 
k 



(8) 



with 



rr d r,,^^ / I d f d 

Cdis = -VcTT + K{E) _ + -_ r— 
oz \ oz^ r or \ or 



One needs to solve a complete triangular-like set of 
coupled equations since a given nucleus can only be 



a lighter. Quantities in this equation are functions 
of spatial coordinates (not time, steady-state being 
assumed) and of kinetic energy per nucleon (energy 
for short) since this is the appropriate parameter to 
be used, as it is conserved in spallation reactions. 



j3 disintegrali 




II mugiielic iiilioiiiogeneities 



Acceleration by shock waves 



Figure 2: Schematic view of our Galaxy as well as 
all propagation steps included in our model. 
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2.3.7 Solar modulation 

The cosmic rays that we detect on Earth had to pen- 
etrate the Solar cavity, a process by which they lose 
energy. This phenomenon is called Solar modulation 
(see for a review), and may be pictured as fol- 
lows. The Sun emits low energy particles in the form 
of a fully ionized plasma having v ~ 400 km s~^. 
This so-called Solar wind shields the Solar cavity 
from penetration of low energy GCR. It was first 
studied by Parker who established the evolution 
of the flux in the Solar cavity. Analog to the prop- 
agation in the Galaxy, it is a diffusion equation in 
a quite different geometry (spherical). For practical 
purposes, this equation can be solved numerically, or 
one can use the force-field approximation. 

Force-field Perko provided a useful and com- 
pact approximation to the full modulation equation. 
The final result for a nucleus with charge Z and 
atomic number ^ is a mere shift in the total energy 



E 



iTOA , 



/A = E'^/A - \Z\ (t>/A 



(9) 



Here E^'~'^ and E^^ correspond to the top-of- 
atmosphere (modulated) and interstellar total en- 
ergy, respectively. The Solar modulation parameter 
(j) has the dimension of a rigidity (or an electric po- 
tential), and its value varies according to the 11-years 
Solar cycle, being greater for a period of maximal So- 
lar activity. Another equivalent quantity often used 
is $ = \Z\(I)/A ~ ^(j). Once the momenta at the 
Earth p^OA ^y^q boundaries of the heliosphere 

p^^ are determined, the interstellar flux of the con- 
sidered nucleus is related to the TOA flux according 
to the simple rule 



TOA 



(i;TOA) 



$IS (_EIS) 



^TOA 



P 



IS 



(10) 



The determination of the modulation parameter 
$ is totally phenomenological and suffers from some 
uncertainties. As explained in the following, we will 
deal with data taken around periods of minimal Solar 
activity, for which we fixed $ = 250 MV (or equiva- 
lently cj) = 500 MV). 

The effect of Solar modulation may be decoupled 
from the problem of interstellar propagation. Would 
a more careful treatment of Solar modulation be 
needed, e.g. j34|, an interstellar flux can easily be 
obtained from force-field modulated fluxes by de- 
modulation (the forcc-fleld approximation we used 
is reversible). This interstellar flux thus obtained 
could then be used as an input for any other pre- 
ferred treatment of Solar modulation. Of course. So- 
lar modulation induces some uncertainty, but this 
question is still debated, and a rigorous treatment of 
this effect is beyond the scope of this paper and more- 
over has not the same impact on antiprotons/protons 



is greater but as one is mostly interested in setting 
upper limits and drawing exclusion plots on exotic 
sources, it is not worth using complex modulation 
schemes. The other nuclei are less sensitive to the 
exact modulation parameter and to the scheme used. 
Thus, our results should not be too sensitive to them. 



2.4 Modelling the effects affecting propaga- 
tion 

We gave above the basic ingredients that enter in 
propagation models. We now come back to several 
aspects for a more detailed discussion, the purpose 
being twofold: first, some ingredients support endur- 
ing developments and for these, it is not always an 
easy task to choose among the variety of parame- 
terizations and approaches available. Thus, taking a 
sort of step by step walk - leading sometimes to kind 
of zoological approach -, wc will sec how problems 
were tackled in the literature in the past, up to the 
point they reached today. Secondly, as modelling the 
galactic environment is a never ending story, because 
the description is richer and richer as one dives into 
fine structure subtleties, we find it advantageous to 
stress the level of detail we wish to include in our 
model. Finally, some choice will be made for the 
sake of tractability of the propagation problem. 



2.4.1 Matter content, ISM and LISM 

The propagation of Galactic Cosmic Rays is influ- 
enced by the distribution of matter in the disc and 
the magnetic structure of the Galaxy. These two in- 
gredients require a careful modelling. As regards the 
flrst point, the matter content is quite well known 
(see for a review) from in situ observations of 
gas thanks to multi- wavelength surveys. As regards 
the second point, it is more difficult to determine be- 
cause of our particular position in the galactic disc. 
Studies of other spirals are particularly useful in that 
respect. For instance, the presence of a diffusive halo 
or a galactic wind is more easily established in these 
other galaxies by radio observations. 

Local interstellar medium (LISM) As will be 
discussed in Sec. 6.3 some radioactive species of cos- 
mic rays are sensitive to the very fine structure of 
the interstellar medium (ISM), in a region of about 
^ 100 pc around the Solar neighborhood, usually re- 
ferred to as the local interstellar medium (LISM). A 
general review of our local environment can be found 
in The LISM is defined as a region of radius 

< 65-250 pc (the bubble) containing extremely hot 
(~ 10^ - 10^ K) and low density (n < 0.005 cm^^) 
gas, surrounded by a dense neutral gas boundary {hy- 
drogen wall) |]37|] . A finer description of this bubble 
tells that the Sun is located in a local fluff with iVni ~ 
0.1 cm~^, T ~ 10"* K and a typical extension of ~ 50 
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modelling to realize that the local bubble is highly 
asymmetric [|6[ and that several cloudlets are 
present in the bubble (for a schematic representa- 
tion, see Fig. ||). Various models have attempted to 
explain the formation of this local bubble 
but this subject is far beyond our concern. 

Several approximation levels may be considered to 
model this bubble. The first level is that for most 
species of cosmic rays it may be ignored, as the pres- 
ence of the bubble has almost no effect on their prop- 
agation. A next level is necessary for some radioac- 
tive species which are very sensitive to the local envi- 
ronment. In the thin disc approximation, this bubble 
may then be represented by a lower density in some 
disc area surrounding the Sun. In our model, it is a 
circular hole of radius rhoie (to be determined) with 
a null density. It will be called "the hole" for short. 
A third level of approximation is the one proposed 
in , who use a three-layer model. 




Figure 3: Schematic representation (not drawn to 
scale) between the local bubble and the neighboring 
Loop I superbubblc (from pO||). 



Interstellar medium (ISM) Apart from the ra- 
dioactive species, cosmic rays travel for hundreds of 
kpc before reaching us, so that they are not sensitive 
to the kind of structures mentioned above but only 
to the smoothed out properties of the ISM. From 
this coarse point of view, the ISM is a rather homo- 
geneous mixture of neutral hydrogen, helium, molec- 
ular hydrogen, ionized gas and dust (see e.g. p7|) 
located in a narrow disc. The total density of matter 
varies with the distance r, but for the sake of an- 
alytical tractability, we consider a constant density 
of ^ 1 particle cm~'^. This is not a crucial assump- 
tion, as the cosmic rays detected on Earth have only 
probed a region of a few kpc of extension, so that 
this effective number can be considered as the mean 
density in this region. 

To end with the ISM question, a last point deserves 
attention. Gas forms clouds, and stars form with 
gas, especially where dense clouds are. In [Q, au- 
thors have studied the consequence of having most of 



question of a larger production of secondaries in these 
clouds, they could have different magnetic fields and 
then be more or less transparent to cosmic rays. As 
a matter of fact, we would be left with two more 
adjustable parameters (filling factor, transparency), 
and in a certain way, other usual propagation param- 
eters fit to observations would have the meaning of 
averaged effective quantities, analogously to a uni- 
form distribution of source and gas. Anyway, with 
our model, we do not intend precisely to provide bet- 
ter than effective parameters, useful and applicable 
for all species that seem to have the same propaga- 
tion history. This notion of effectiveness and domain 
of validity of our model will often come back all the 
paper long. 



2.4.2 Diffusive halo, regular and erratic magnetic field 

Many aspects of cosmic rays propagation call for 
the existence of an extended magnetic halo (see 
e.g. iQ). First, clear evidences were obtained from 
non-thermal radio emission in NGC 4631 |Q. Two 
components of the galactic magnetic field coexist: a 
regular one (average value about a few /iG, parallel to 
the galactic plane, responsible for confinement) and 
a stochastic one which is responsible for charged nu- 
clei diffusion (as well as diffusive reacceleration), that 
has about the same strength. Tab. |l|, taken from psj , 
summarizes several observational evidences. In the 
last ten years, several studies tried to constrain the 
halo size using a combination of stable plus ra- 
dioactive nuclei ratio measurements. They found 
L - 4 kpc [|o|, i < 4 kpc @, L ~ 4 - 12 kpc [|3 
and L ^ 2 — 4 kpc |^|. As all the authors used 
quite different diffusion parameters, with or without 
wind, reacceleration. and as they took quite differ- 
ent diffusion slopes 5, the trend seems to favor small 
halos. However, the radial 7-ray distribution in our 
Galaxy is rather Hat, which points towards a large 
halo (see e.g. |4^]). To face the problem, |Q have 
recently proposed a more complex diffusion model, 
where some quantities (galactic wind) have a spatial 
dependence. We feel that such an approach becomes 
unavoidable with the present development of cosmic 
rays studies. However, we wish to mention two ar- 
guments against the claim of inconsistency advanced 
by these authors to motivate their model: first, we 
showed in |4^, |4| that considering B/C ratio in diffu- 
sion/convcction/reacceleration models, all halo sizes 
are possible. Second, focusing on radioactive species, 
we showed that if we agreed with usual results 
for a homogeneous LISM, the existence of a local 
underdensity leaves the choice of L almost free (too 
large halo sizes, i.e. L > 12 kpc are excluded). Thus, 
from our point of view, the discussion about the halo 
size is far from being closed. 

Finally, it must be kept in mind that halos are dy- 
namical objects. In particular, there seems to be a 
transition zone (thick disc) between the more quiet 
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VPS ['?| 


> 3 kpc 
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Gradient in CR density 








a- data from SAS-2 
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~ 3 kpc 




b- data from COS-B 
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a,- p + p 
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> 1 kpc 




b- Inverse Compton 


yes 


- 10 kpc 



Table 1: Estimations of the cosmic ray halo extension (adapted from |{45|| ). 



from the clouds, hot ionized gas and stars kinemat- 
ics. The halo could be substructurcd in chimneys, 
galactic fountains and galactic holes [Q. The im- 
portance of this substructure for cosmic ray studies 
is an open question. 



2.4.3 Standard sources in the disc 

One easily imagines that the location of sources plays 
an important role for the density of all species, pri- 
maries, secondaries and radioactive nuclei. More- 
over, it is intimately related to the question of valid- 
ity of the stationary hypothesis in cosmic ray prop- 
agation models. Nearby sources are for example in- 
voked in to explain the observed behavior at PeV 
energy, i.e. the knee. If one goes one step further into 
this line of thought, nothing would guarantee that 
even at low energy the steady state approximation 
is justified. It could be that propagation parame- 
ters provided by stationary models are just effective 
parameters, that would need to be reinterpreted in 
terms of time dependent propagation models. Such a 
view could seemingly solve some present problems of 
cosmic rays (Maurin et ai, in preparation). Anyway, 
this paper focuses on stationary models. 



Radial distribution of sources q{r) Measure- 
ments of galactic 7 rays in the seventies have raised 
the question of the radial distribution of cosmic rays. 
This distribution is needed in order to evaluate the 
resulting 7 emissivity at different galactocentric lo- 
cations. The first distribution used was that of Ko- 
daira | |5^ following the radial distribution of super- 
novae which is also close to that of pulsars. This 
is consistent with the present picture of cosmic rays 
where supernova: provide the energetic budget and 
mechanism to accelerate nuclei. The description of 
the galactocentric distribution has been steadily im- 
proved thanks to new observations of pulsars and 
supernovae. We take here the distribution of Case & 
Bhattacharya [p5[ w hich is an improvement of their 
earlier analysis [p6| and happens to be closer to the 



All these distributions could be compared to flat dis- 
tribution, i.e. q{r) — 1, that is also widely used 
in one-dimensional propagation models These 
three radial distributions are 

• Model a: flat distribution 

q{r) = 1 . 

• Model b: Case & Bhattacharya 



• Model c: Strong & Moskalenko p7 



q{r) 



r 

8^ 



0.5 



cxp 



-1. X 



(r- 8.5) 



They are displayed in Fig. || along with the metallic- 
ity gradient discussed below. 

We emphasized in that the difference between 
a flat and a more realistic q{r) was a mere rescal- 
ing of the propagation parameters. However, this 
difference is not so important if one considers these 
parameters as effective and uses them in a consistent 
framework. To be more specific, one should make 
the same hypothesis on q(r) when i) determining 
the propagation parameters from the observation of 
some cosmic species and ii) predicting the abundance 
of some species from these parameters. On the other 
hand; if one is interested in a physical interpretation 
of these parameters and wants to deduce some galac- 
tic property or to predict some quantity which docs 
not depend only on the diffusion properties (such 
as proton- induced 7-ray production), one should be 
careful to consider the right radial distribution. 



Metallicity gradient It is also expected that the 
composition of SN should depend on its position in 
the Galaxy. Indeed, the inner Galaxy is richer in 
heavy stellar material than the outer Galaxy; this 
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Model b 
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Metallicity gradient 
— not Included (ttiln lines) 
uded (ttiick lines) 
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r [kpc] 

Figure 4: Radial distribution of the sources. The 
thin hues correspond to the p and He component, 
with the three types (a), (b) and (c) discussed above. 
The thick Hnes correspond to the other species, for 
which the distribution id modified by the metaUicity 
gradient, according to Eq. (p^). 

showed a gradient for O /H from observations of ion- 
ized nebulae in galaxies like M33, M51, and MlOl, 
but later works observed this trend in our Galaxy 
for many other abundances (sec for a review and 
Sec. 2 of ll^). Several recent observations (see Tab. 1 
of Q for a compilation and for latest results) 
lead to very similar conclusions for the metallicity 
gradient. It can roughly be parameterized, for all 
ion species X but He, as 



d[X/H] 
dr 



-0.05 dcx kpc 



(11) 



where [X/H] is defined as Zo,9io(X/H) — Zo,9io(X/H)q. 
Associated with the radial distribution, this gives an 
additional factor 



little effect on the origin of sources that contribute 
to the flux detected at Earth. However, we are left 
once more with the question of the meaning of the 
propagation parameters we will derive. 

Nevertheless, if a spiral structure is retained, the 
other ingredients should follow as well (magnetic field 
structure, gas distribution,. . . ); such a modelling is 
beyond the scope of this study, though this is prob- 
ably an important point. Joint progress in several 
research fields (diffusion mechanisms, better knowl- 
edge of galactic structure, etc) associated with nu- 
merical developments will certainly allow to cope 
with these questions in a couple of years. Anyway, 
we will not discuss further this point in this paper. 

Source spectrum The energy spectrum of the 
particles emitted by the sources is determined by the 
acceleration process at work. There is a strong belief, 
based on theoretical work, that this energy spectrum 
dQ/dp is a power-law in rigidity |^^. For a species j, 
the differential spectrum in energy dQ^ / dE = Q\TZ) 
is then given by 



n 



(3 \1GY 



(13) 



in our Eq. (||) (see details in ) , where the value of 
a is still debated (grossly from about 1.5 to 2.5 jl^). 
We will refer to this case as the pure power law spec- 
trum. Among the various acceleration models, some 
of them are able to produce various spectra for dif- 
ferent nuclei JSTf , but also a break in slopes around 
PeV energies (see references in [Q). 

Because of collective effects, one can also consider 
a possible deviation from this power-law at low en- 
ergy, and we choose 



n 



1 GV 



(14) 



We refer to this as the modified spectrum. 



qir) iQ-O-OHr-s.) X q(r) , 
except for H and He. 



(12) 



Spiral structure Because of the presence of spi- 
ral arms, there is also an angular dependence of the 
source distribution. As showed, SNe II and lb 
are likely to follow the gas distribution. As also con- 
firmed in , SNe are seemingly correlated with the 
spiral arms position. The structure can be found in a 
recent meta-analysis of the Milky Way where the 
structure is likely to resemble a 4 logarithmic spiral 
arms (see its Fig. 2). Notice that, despite their pos- 
sible usefulness, such models already fall to describe 
the arms outside the range 3—13 kpc. Anyway, it 



2.4.4 Exotic sources in the halo 

Cosmic rays may also be created by the pair anni- 
hilation of exotic particles - such as the ncutralino 
predicted by supersymmetry (SUSY) -, as discussed 
in Sec. 7.3, or by the evaporation of primordial black 
holes - hereafter PBH - as discussed in Sec. |7.4 Both 
populations are supposed either to fill (SUSY) or to 
follow (PBH) the dark matter profile, which can be 
assumed to be a spherical isothermal profile, i.e. for 
each species pi{r, z) = pf x f{r, z) with 



fir,z) 



(15) 



For definiteness, the dark matter halo has been as- 
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uncertainties on Rc and the consequences of a possi- 
ble flatness have been shown to be irrelevant at least 
in the PBH case [lei. 



The source density in the Solar neighborhood is 
about p® = 0.4 GeV cm~^. For PBH the normahza- 
tion /OpBjj is very different (far less) and can be con- 
strained through CR signatures (see Sec. 7A). This 
is at variance with the SUSY case where the SUSY 
space parameter is constrained instead of source den- 
sity (see Sec. 7.3). 



Effective source term Finally, the production 
rate, i.e. the effective source term is different in the 
two cases: it is proportional to p%g[j x /(r, z) for 
cosmic rays evaporating from PBH, but to (p® x 
f{r,z))'^ for those coming from neutralino annihi- 
lations. The resulting source profiles arc shown in 

(kpc) 



Dark Matter 
Isothermal profile 
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2h=0.2kpc 
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. (ISM) 
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Convection Vc 
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( iixial symmetry around zl 



Figure 5: Schematic view of our Galaxy: exotic pri- 
mary CR production follows either the dark matter 
profile (PBH, solid line) or the square of the dark 
matter profile (SUSY, dashed line). They have been 
arbitrary rescaled. 

Fig. ^, compared to the typical radial extension of 
the Galaxy. One easily imagines that sources closer 
to the galactic center will contribute more to the de- 
tected flux, even more in the SUSY case than in the 



PBH case (see also discussion in Sec. 7.2 ) 



2.4.5 Catastrophic losses; destruction on interstellar 
matter 



The overview in Sec. 2.3.2 emphasized the role of nu- 



clear reactions and decay in the propagation of cos- 
mic rays. The importance of these events is not the 
same for all species, and one should compare the typ- 
ical propagation time ~ 20 Myr to the effective time 
of the process involved. The most important loss is 
from spallations, i.e. inelastic interaction of the nu- 
cleus with the matter content of the thin gaseous disc 
(as a crude estimation, H : He : CNO : Fe « 10 : 1 : 
0.1 : 0.01, about the same proportions as for cosmic 



of secondary species: LiBcB and sub-Fe nuclei (Sc, 
Ti, V) are absent in sources but they are observed 
with the proportion B/C w sub-Fe/Fe w 0.1 in cos- 
mic ray radiations; they were formed by spallations. 
Of course, the heavier nuclei have larger cross sec- 
tions and break more easily during their wandering 
through the Galaxy. This preferential destruction in 
flight is, as a matter of fact, responsible for part of 
the nuclear enrichment in heavies observed till PeV 
energies (p3, see also Sec. 6.4). 



When a cosmic ray nucleus A of energy Ea inter- 
acts with a nucleus H or He (for simplicity, we write 
the formulae for H only) at rest in the interstellar 
medium, a nuclear reactions may occur, producing 
some daughter species B, C, ... of different energies 
Eb, Ec, ■ ■ ■ Actually, it turns out that the kinetic en- 
ergy per nucleus E is approximately conserved dur- 
ing the reactions, so that all the information we need 
is contained in the set of differential cross sections 
d(7A+H->B+c+.../dE. From these, the total inelastic 
cross section for the species A may be obtained by 
summing over the final states B, C, ... The produc- 
tion cross section for the species B can also be found 
by summing over the initial states A. 

On a purely theoretical point of view, the compu- 
tation of the duA+H^B+c+.-./dE is a well-defined 
but hard to solve nuclear physics problem (see 
for a review). A nucleus- nucleus collision is generally 
an out of equilibrium process that need to be dy- 
namically treated. At a few hundreds of MeV/nuc, 
the system expands and the equation of state of nu- 
clear matter enters in the spinodal region (region 
of great mechanical instability) : small perturbations 
are amplified until they reach the size of the sys- 
tem. At GeV energies, one is left with several nu- 
clei (multi- fragmentation) ; a statistical description 
can be used (the number of fragments does not de- 
pend on energy but only on the nucleus involved). 
More specific processes, such as single and double 
nucleon removal are indeed the dominant processes 
for many nuclei at these energies. Purely hadronic 
processes arise at a few hundreds of GeV; more 
and more resonances are excited until quark frag- 
mentation becomes dominant. Among microscopi- 
cal approaches ||7^, 0, one could quote probabilist 
models (statistical multi- fragmentation) , dynami- 
cal models (quantum molecular dynamics) and ki- 
netic models (Boltzmann/Vlasov-Uehling-Uhlcnbeck 
equations). However, none of these models enables 
a quick and precise evaluation of the desired cross 
section. At most, they can reproduce data with only 
a reasonable agreement. This approach is certainly 
promising as empirical formulae have now probably 
reached their highest point of refinement and devel- 
opment. 

From the experimental point of view, the 
accurate determination of the cross sections 
duA+H^B+c+... / dE is difficult, as it requires the 
identification of all the nuclei produced in the re- 
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A is much easier to measure, as one has only to know 
that the reaction actually occurred. This is the rea- 
son why the question of total reaction cross section 
is usually separated from the problem of production 
cross section of secondaries. They will accordingly 
be considered separately in the next two parts. 



Semi-empirical formulae |72| The first ingre- 
dient is the reaction (or inelastic) cross section 
f^inciiE), giving the probability that a cosmic ray nu- 
cleus of energy E undergoes a nuclear reaction with 
an interstellar nucleus. It should ideally be given by 
a fundamental approach, but as discussed above, this 
program has not been fulfilled yet. Simpler consid- 
erations lead to semi-empirical formulae which turn 
to be quite powerful in reproducing all existing data 
on nuclear reactions. In the simplest approach, the 
cross section is proportional to the geometrical area 
of the nucleus (Tri?^^^) that scales as ^ A^l'^ . In 
1950, Bradt & Peters have proposed a first correc- 
tion to take into account the overlap of the two nuclei 
wavefunctions, 

CTinci 7rro(At^,g„t + ^proj ^ Oo) , 

where tq is the nuclcon radius and the overlap 
parameter (or transparency). 

In the eighties, Letaw and coworkers provided a 
refinement which is still widely used, in particular 
for p -I- A reactions. At high energy, the cross section 
tends to the asymptotic energy independent form 

a™i(mb) = 45^°-^[l + 0.016sin(5.3 - 2.63 In^)] 

At lower energy, the behavior becomes energy depen- 
dent: one finds a slight hollow around 200 MeV/nuc, 
followed by a sharp enhancement and a sharp de- 
crease till ^ 20MeV/nuc; 



HE 
incl 



HE 
incl 



(mb) 



X [1 - 0.62 exp(-i;fc/200) sin(10.9£;^°-2^)] . 



In the above formula, denotes the kinetic energy 
per nuclcon in MeV/nuc. More generally, there are 
some necessary corrections for light nuclei (due to 
their particular nuclear structure) and the reactions 
He+A are evaluated through an energy dependent 
rescaling factor [crHe/o'p](£'fc). Later on, in 1993, Si- 
hver and coworkers, gave a formulation closer to the 
geometrical picture given above, to take also into ac- 
count A-l-A reactions, thanks to new experimental 
data. Finally, in 1996, Wellish and Axen introduced 
refined functions in Sihver et al. formulae in order to 
be able to describe the whole data set with a 2% ac- 
curacy from 6.8 MeV/nuc to 9 GeV/nuc (more than 
1400 data points, but they only retain points whose 
accuracy is better than 4%). 

Universal parameterization [ItsII As was shown 



fit to most measurements of nucleus-nucleus inelastic 
cross sections is currently the universal parameteri- 
zation proposed by p3|. It is given by 



r2 (a^/^ 

'^0 l^^proj 



A 



1/3 

target + "E 

B 



(16) 



1 - Rc 



E, 



Here, rp = 1 fm and Ecm denotes the kinetic energy 
per nucleus in the rest mass frame. Rc is a parameter 
introduced for light systems in order to keep the same 
formalism for all reactions. The term 



1/3 



E, 



describes several effects: 



• S 



AtAp 



A^^^) is the asymmetri- 



cal mass term that is related to the overlapping 
volume of the colliding system; 

• Ce = D[l-c^x>{~E/Ti)] - 
0.292 exp (-£;/792) cos(0.229£;°-453) models 
the energy dependence at intermediate and 
large energies, mainly through transparency 
and Pauli blocking {E is the colliding kinetic 
energy per nucleus in MeV/nuc). The number 
D is related to the system density and as Ti, it 
is adjusted once for all for almost all reactions 
with a more careful treatment and specific 
adjustment for light systems [A <A). 

• Finally, the last term in 5e takes into account 
isotopic dependence of cross sections. 

The remaining terms in Eq. ([l6|) , respectively 
lAAZpZt 



B 



R 



and 



1.2 U 



1/3 



A 



1/3 



R 



E, 



1/3 



correspond to the coulomb barrier (that depends on 
the energy) and the radius to evaluate the barrier 
height. In the latter term, is the equivalent radius 
of the hard sphere and is related to rrms,i thanks 
to Ti = 1.29 rrms,i (this last quantity is obtained 
directly through experiment ||75|| ). 



Experimental data and accuracy [76] As an il- 
lustration. Fig. ^ shows these inelastic cross sections 
for the reactions (p,IIe)-|-^Be, (p,He)+^^C along with 
the experimental data points. We see that this pa- 
rameterization reproduces well the non-trivial fea- 
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Figure 6: Comparison of total inelastic eross sections from universal parameterization to data |7^: empty 
triangles are from Bauhoff (1986) and Carlson (1996), full ones from Tanihata et al. (1985), and diamond from 
Bobchenko (1979). 



stressed that no new measurement has been made 
since 1985 for the reaction of protons on light nu- 
clei, and that many isotopes have even never been 
measured to date. Moreover, the helium induced re- 
actions are even more difficult to probe experimen- 
tally. Some of the results concerning the abundance 
of some species of cosmic rays rely on the faith that 
the parameterization is valid for these nuclei. Con- 
sidering that when data exist, the parameterization 
gives an accuracy of a few %, it is reasonable to con- 
sider that the above parameterization is 2 — 5% ac- 
curate for proton induced reactions and 10 — 20% for 
He induced reactions. 



2.4.6 Secondary production from spallations on the 
interstellar matter 

When a nucleus undergoes a spallation, it may gen- 
erate a large variety of lighter nuclei, which are not 
all present in standard sources. These nuclei may 
typically be classified as ^H, ^He (p and He in- 
duced), LiBeB (CNO induced) and sub-Fe (Fe in- 
duced) . They give important clues to understand the 
characteristics of galactic propagation. As another 
example of the importance to have a good knowl- 
edge of these production cross sections, one can also 
mention cosmogenic nuclei (due to the penetration of 
GeV CR protons in extraterrestrial bodies) which al- 
low an estimation of the stability of cosmic ray fluxes 
in the past billion years |7^] . To study the nuclei pop- 
ulation N2 and to draw an inference for propagation 
models, we need to now, as accurately as possible, 
processes such as 



In a more rigorous way, one should also take into ac- 
count contributions from heavier species of the ISM 
(e.g. A^i-|-CNO). However, these cross sections are 
only poorly estimated, not to say not known; more- 
over, these have been correctly included for antipro- 
ton production but were shown to be negligible con- 
tributors. 



Semi-empirical and empirical formulae [78| In 



the literature, there exists basically two fast proce- 
dures to evaluate the spallation cross sections cor- 
responding to reactions ([T^). They were provided 
by cosmic ray physicists because the energy range 
in which these reactions are needed is generally not 
of great interest for nuclear physicists. Silberberg 
& Tsao approach has well-motivated theoretical 
grounds and takes advantage of some observed regu- 
larities in i) the mass difference between target and 
daughter nuclei (AA = Aargot - ^produced); ii) the 
ratio between the number of neutrons and protons 
in the daughter nucleus. 

In the nineties, many new data issued from various 
targets were obtained, and this semi-empirical pa- 
rameterization gave only a ±35% accuracy for these 
new data. It led Webber and coworkers to de- 
velop a new approach fully based upon experimen- 
tal regularities. In this parameterization, having 
remarked that production of different isotopes has 
similar E dependence, parameters were fit to one ki- 
netic energy per nucleon, and then extended to other 
energies, 

af+^-(Z„^, ;Z/,^^;£;) =ao(Z/,Z,) 
xh{Zf,Af,Z,,A,)h{E,Zf,Z,) . 
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ferred to quoted papers), but only briefly comment 
the three terms that enters in this equation: the 
first one describes the dependence on the fragment 
charge (the production is exponentiaUy suppressed 
far from the stabiUty valley) and contains in partic- 
ular two parameters {azf and Azf) that are chosen 
to fit the data. Notice that these parameters, thanks 
to new measurements, were updated in a later ver- 
sion of their code; we implemented in our propa- 
gation code parameters as given in Tab. V of | pl[ . 
The second term describes isotopic distribution of 
fragments for a given species. This is independent 
of energy as long as we are not too far from the 
stability valley. The third term is the energy de- 
pendence and is only related to charges involved. 
Added to these three terms, the particular single and 
double nucleon-removal reactions are also taken into 
account, and He- induced spallations arc tabulated 
from 1^ formulae. 

Some remarks about these parameterizations 

Actually, for our purpose, the most suitable set of 
formulae is given by Webber and coworkers' empir- 
ical formulae. It seems that the more general ap- 
proach of Silberberg and coworkers is more useful to 
propagate Z > 30 heavy nuclei (see also Q) as the 
corresponding cross section have not been measured 
yet. It was found in |Q that for the existing data 
points, the dispersion obtained with Webber's code is 
better than the dispersion using theirs. However, the 
two approaches have probably reached their latest 
development. As remarked recently in plf , any re- 
finement of these formulae asks for complicate mod- 
ifications of the present formulae. If one takes into 
account isospin effect, odd-even effect, neutron and 
proton stripping due to coulomb barrier, it would re- 
quire an adjustment nucleus by nucleus. Webber's 
code is only valid for Z > 3: light systems require 
further modelling because they behave very differ- 
ently from other nucleus (see e.g. Q). To improve 
these codes, some authors have focused on peculiar 
single and double nucleon removal (coulomb frag- 
mentation) that are dominant for most secondary 
productions jsH) . A further point of discussion is the 
reliability of the straight-ahead approximation, i.e. 

/>oo 

/ (3'N^{T')cj^^{T,T')dT' l3N''{T)a''^{T) . 
Jo 

In practice, this corresponds to the fact that each 
fragment has the same kinetic energy per nucleon 
than has the parent nucleus. If this approximation 
seems quite justified at high energy (above a few 
GeV/nuc), its relaxation may affect, for example, the 
B/C ratio around 1 GeV/nuc |Q. In our propaga- 
tion code, we use Webber et al.'s code that is not only 
very well suited, but also widely used by almost all 
cosmic ray physicists (for the nuclei we deal with). 

Experimental data and accuracy [^t!] At the 



the reactions involving cosmic ray nuclei have been 
measured at least at one energy. As an illustration, 
we give in Tab. ^ the cross sections for the reac- 
tions involving ^^C. It illustrates the importance of 
single nucleon removal (^^C ^^C, ^^B) as men- 
tioned above. Double nucleon removal (^^C ^"B) 
is less marked but visible. Let us notice also a maxi- 
mum around ^Li that corresponds to half of the ini- 
tial nucleus. This table also emphasizes the multi- 
fragmentation character of spallations. 



IXcdCLlUli 


Cspal ± AfTspal 




fmb) 


C 


29.2 ± 2.5 


lOr< 


3.6 ± 0.5 


9p 


n 9/1 -1- n n ^ 




0.12 ±0.05 


"B 


27.7±0.7 


10b 


12.3 ±3.0 


«B 


0.44 ±0.04 


i"Bc 


4.2 ±0.6 


9Be 


6.7 ±0.9 


^Bc 


10.1 ± 1.3 


*JLi 


0.25 ±0.06 


8Li 


1.47 ±0.23 




12.5 ± 1.8 


6Li 


19.8 ±2.7 


«He 


0.87 ±0.31 


^He 


159. ±21. 


^Hc 


24.8 ±3.2 


^H 


88. ±31. 


2H 


138. ±41. 


iH 


143. ±42. 



Tabic 2: Cross sections of the fragmentation of ^^C 
on a proton target at 3.66 Gev/nuc (adapted from 
Korejwo et al, 1999) 0. 



Tab. ^illustrates another aspect of spallations in cos- 
mic ray propagation: the output is from Webber's 
spallation code and has been weighted by the rela- 
tive abundances of sources to isolate the most im- 
portant contributions. This table underlines the fact 
that, if most secondaries are created through spal- 
lation of the most abundant species, all the others 
come from small contributions of a great numbers 
of nuclei. Finally, the present data accuracy ranges 
from ^ 2 — 3% to ^ 20% depending on the pro- 
duced fragment. The old activation measurements 
seem to be unusable, and despite numerous experi- 
ments since the beginning of the eighties, not all the 
nucleus have been at least measured twice. Recently, 
Flesh et al. (1999) js^ observed systematic devia- 
tions from Webber measurements. It is thus not an 
easy task to estimate the accuracy of the above pa- 
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Total 


94.4% 


92.1% 


96.5% 


93.0% 


93.7% 



Tabic 3: Fraction in % of the isotopes contributing 
to the formation of secondary B and Be, (when this 
fraction is greater than 1 %), at 1.8 GeV/nuc. The 
contribution of each channel is weighted by the prop- 
agated abundance of the father nucleus. 



particularly important in our studies, one can esti- 
mate the accuracy to lie in the range of about 5—10% 
for p induced reactions and 10 — 20% for He induced 
reactions. It is worth noting that production cross 
sections will probably be the limiting factor of cos- 
mic ray studies in the near future, when high statistic 
experiments (pamela, AMS) will obtain new very ac- 
curate data on CR fluxes. This is already the case 
for antiproton predictions (see jssj and Sec. |6.l| ). For 
secondary nuclei such as B, sub-Fe, consequences of 
these errors are difficult to estimated, because of the 
complex reaction chain involved (see |^ for a dis- 
cussion). 



Ghost nuclei There remains a last point which is 
only rarely addressed. As mentioned above, a com- 
plete grid of nuclear reactions is needed to imple- 
ment the cross sections, including short-lived species 
which are not seen directly but that contribute to 
some other species through their decay products. As 
an example, let us consider the production of ^Be. 
One has to add all the contributions from all heavier 
nuclei, including the specific contribution from inter- 
mediate reactions such as 



(p,He) 



^Li 
"Li 



(ti/2 = 178 ms) 
fi^+ 2n (ti/2=8.6 ms) 



^Bc (Br = 0.49) . 
^Be {Br = 0.041) 



The nuclei ^Li and ^^Li are called ghost nuclei: they 
contribute indirectly to the flux, but do not have to 
be propagated as their half-lives are negligible in the 
face of propagation time scale. This will be the case 



that there is a gap in the lifetimes of nuclei between 
10 kyr and a few Myr, so that all nuclei having ^ 
1 kyr could be considered as ghosts. Among them, 
those which have too short a lifetime are not seen in 
the experiments and are already accounted for in the 
measured cross sections. We emphasize that, would 
we be interested in supernovae explosions, the nuclei 
relevant would not be the same because of a much 
shorter typical evolution time. 

In conclusion, even if it will be transparent in the 
rest of the paper - and also in almost all propagation 
papers -, one has to keep in mind that each time we 
evaluate production cross section, in the true com- 
putation intervenes the quantity {Br(X — > j) is the 
branching ratio of X into j) 



^effective 



^direct 



ghosts X 



Br{X ^ j) 



Nuclear grids for all species (including very heavy nu- 
clei) can be found in j90|. Thanks to recent nuclear 
compilations |Q, we reevaluated this data grid (and 
the associated ghost nuclei) for Z < 30; we found 
many new ghost nuclei or modified branching ratio. 
However, these new nuclei are generally far from the 
stability valley, such that their production cross sec- 
tion is clearly very small. We checked that for most 
species, this completion have a negligible effect, so 
that the [@'s table of ghost nuclei is sufficient and 
does not need to be further revised R- 



Production of antiprotons by spallations 

Among the nuclei that are created by spallations, 
the antiprotons require a specific treatment, as their 
creation relies on the conversion of collision energy 
into mass and not on the mere breaking of a heav- 
ier nucleus. In order to evaluate the contribution 
from the p — Hism interactions - the dominant pro- 
cess at first approximation - the Tan & Ng param- 
eterization based on experimental data is used p^ . 
Collisions that involve heavier nuclei are more dif- 
ficult to obtain. The contribution p — Hcism has 
been considered firstly by means of a simple geomet- 
rical approach However, using a more sophisti- 
cated nuclear Monte Carlo treatment, it was noticed 
in j9^ that this reaction does not only enhance the 
antiproton flux as a whole but also change its low en- 
ergy spectrum, mostly for kinematical reasons. By 
the way, these authors also considered He — pisM, 
He — HeisM as well as CNO contributions, the lat- 
ter being negligible (see their Fig. 7) . Unfortunately, 
very few experimental data are available on antipro- 
ton production cross sections in nuclear collisions 
and a model-based evaluation - such as the dtunuc 
program |9^ - seems to be unavoidable. This is what 
we will use in the following. 



■^A comnlete table of Erhost nuclei for Z < .SO can be found 
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2.4.7 Antideuteron production 



When an antiproton and an antineutron are formed 
by a reaction, they may merge into an antideuteron 
nucleus. The calculation of the probability for the 
formation of an antideuteron proceeds in two steps. 
We first need to estimate the probability for the cre- 
ation of an antiproton-antincutron pair. Then, there 
is some probability that those antinucleons merge to- 
gether to yield an antinucleus of deuterium. 

As regards the first step, the differential probabil- 
ity for the production of a single antiproton or an- 
tineutron (resp. Vp or Vn) is known for each initial 
process - spallation, neutralino annihilation or black 
hole evaporation - and a first guess would be that 
the production of two antinucleons is proportional to 
their product. This so-called factorization hypothe- 
sis is fairly well established at high energies. For spal- 
lation reactions, however, the bulk of the antiproton 
production takes place for an energy ^/s ~ 10 GeV 
which turns out to be of the same order of magnitude 
as the antideuteron mass. Pure factorization should 
break in that case as a result of energy conservation 
and needs to be shghtly adjusted. We have therefore 
assumed that the center of mass energy available for 
the production of the second antinucleon is reduced 
by twice the energy carried away by the first antin- 
ucleon. This yields the following factorization of the 
probability to form a p-n pair (sec ||95|| for a more 
detailed discussion) 



1 

2 



VpiV^, kp) VniV^ - 2Ep, kn) 



Once the antiproton and the antineutron are 
formed, there is a finite probability that they com- 
bine together to give an antideuteron. The coales- 
cence function C(kp — kn) describes the probability 
for ap — n pair to yield by fusion an antideuteron, as 
a function of the difference of the initial momenta. 
An energy of ~ 3.7 GeV is required to form by spal- 
lation an antideuteron whereas the binding energy of 
the latter is B ^ 2.2 MeV. The coalescence function 
is therefore strongly peaked around kp — kn =0. 

We considered that the antinucleons merge to- 
gether if the momentum of the corresponding two- 
body reduced system is less than some critical value 
-Pcoai- That coalescence momentum is the only 
free parameter of this factorization and coalescence 
scheme. As shown in |9^, the resulting antideuteron 
production cross section in proton-proton collisions 
is well fitted by this simple one-parameter model. 
The result, for a typical value po = 160 MeV as 
well as the comparison with the antiproton yield is 
displayed in Fig. 

New calculations for such astrophysical appli- 




antideuteron energy (GeV) 

Antideuterons spectrum 

Figure 7: Upper plot: antiproton differential spec- 
trum obtained with 1.9 x 10^ u quark jets generated 
by PYTHIA at 100, 75, 50, 25 GeV. Lower plot: an- 
tideuteron spectrum obtained in the same conditions 
with a coalescence momentum po = 160 MeV. 



progress in the diagrammatic approach to the coa- 
lescence model [P6|. 



2.4.8 P and EC decay 

A second catastrophic loss to add to inelastic reac- 
tions is the decay of unstable nuclei. One can distin- 
guish two decay processes: /3 decay which is spon- 
taneous and electronic capture decay (hereafter de- 
noted EC for short) which can happen only if an 
electron of the interstellar medium has been first at- 
tached in the K shell (so that the presence proba- 
bility of the electron is finite in the region of the 
nucleus). 



[3 decay The ubiquitous example of this case is 
i"Bc-^i"B whose rest half-life is 1^/2 = 1.51 Myr. 
The ratio ^°Be/^Be observed in cosmic rays is smaller 
than the ratio of two normal secondaries coming from 
the same primary progenitors, which is merely the ra- 
tio of the production cross sections. This has often 
been interpreted as an evidence for an extended dif- 
fusive halo where the cosmic rays would spend most 
of their time. This interpretation may be mislead- 
ing and could lead to a wrong intuition pO[ | . One of 
the main interest of these nuclei is that the short- 
lived unstable species propagate only very locally, 
and thus are very sensitive to the local environment. 
In particular, they may be used to give an indepen- 
dent evidence for the existence of a loca l undcrdcnse 
bubble, as shown in |5^] (sec Sec. |6.3| for more de- 
tails). 
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agated as pure /3 unstable. Those with half-time 



shorter than ^"'C are ignored (see Sec. 2.4.6 ). In par 



z 


Nucleus 


Daughter 


(error) 


4 

6 

13 

17 

26 


i"Bc 

iAi 

iFe 


26Mg 

?iAr 
iNi 


l.Sl'^y (0.06) 
5.73'^5"' (0.04) 
O.giMyr (0.04) 

Q_3Q7Myr (0.002) 

1.5^^' (0.3) 



Table 4: Pure (3 unstable isotopes (1 kyr < ti/2 < 
100 Myr) from a propagation point of view (see 
for details). 

ticular, we have checked that the EC mode can be 
neglected for ^^Al and ^^Cl but not for ^-^Mn and ^^Ni 
(EC half-life determination for these two last species 
is particularly difficult, see resp. [p^ and [p8[). 



Electronic capture 



99 



A second kind of un- 



stable species is given by nuclei decaying under EC 
process. A nucleus such as ^^Ni is formed during 
last stages of stellar nucleosynthesis; it decays as 
^^Ni— >^^Co in = 80 kyr once it attaches an elec- 
tron in its K-shell. Studies of this species allows for 
example to show that there was a delay of at least 
~ 10^ yr between nucleosynthesis and acceleration. 
Attachment and the converse process, electron strip- 
ping, arc crucial processes that determine the effec- 
tive lifetime of the EC unstable nucleus. Notice that 
because of these processes, the energy dependence of 
lifetimes is more complex in the EC mode than in 
the P mode. 

There are two ways to attach an electron, radia- 
tive and non-radiative process: after attachment of 
a free electron, the momentum and energy balance 
are restored by the emission of a photon or the re- 
coil of the nucleus respectively. Radiative capture 
is dominant around a few hundreds of MeV/nuc in 
hydrogen. Non radiative capture scales as Zf and 
dominates for the heavier species of the ISM. Only 
a few measurements are available for the latter pro- 
cess, and they suggest that this contribution is any- 
way less than a few %; it is neglected in this work. 
It must be said that these processes suffer from the 
same limitations that the production cross sections, 
i.e. they are not studied by atomic physicist, but by 
cosmic ray astrophysicists. The results widely used 
nowadays have been derived in the seventies and to 
our knowledge no updates have been made at GeV 
energies. 

Fig. H displays the rate of attachment and strip- 
ping for several charges. Some general conclusions 
can be drawn from this single result: i) comparing 
^attach to the general propagation time ~ 20 Myr, 
only low energy and large Z nuclei are likely to at- 
tach an electron and subsequently decay; ii) com- 
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Figure 8: Energy dependence of attachment and 

stripping rates Tattach = l/("ISM-f -O-attach) et Tstrlp = 

l/(?T'iSM-w.crstrip) for several charges. 



is captured, decay is much more likely than strip- 
ping; iii) consequence of the previous point is that 
for Z < 30 species, nuclei are always at most singly 
attached. 

Actually, as attachment arises only when there are 
available electrons, one should take into account the 
local properties of the interstellar medium (LISM) in 
all studies involving EC-species (as for radioactive /3- 
decay, for different reasons). These attachment and 
stripping cross sections arc included in the propa- 
gation code used for the work presented here. We 
should mention that these effects may be neglected 
in most cases, e.g. for B/C studies. 



2.4.9 Convection and reacceleration 

The ingredients discussed here were introduced with 
the first theoretical development of cosmic ray 
physics, but they were only taken into account quite 
recently in the propagation models. 



Convective wind It has been recognized for a 
long time that a thin disc configuration would be dis- 
rupted by cosmic ray pressure |10C]. It can be stabi- 
lized by the presence of a halo, but further considera- 
tions imply that this halo would not be static either. 
The stellar activity and the energetic phenomena as- 
sociated to the late stage of stellar evolution may 
push the interstellar plasma and the magnetic field 
associated with it out of the galactic plane. The net 
result would be the presence of a convective current 
directed outwards from the galactic plane and called 
galactic wind, which adds a convective term to the 
diffusion equation. The properties of this wind may 
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hydrodynamical (MHD) [102| methods 



Consequences of a wind have been first investi- 
gated by Ipavitch [103|. Since then, it has been ob- 
served in other galaxies through a flattening of the 
electron induced radio spectrum around 1 GHz in the 
halo 1 104 that could be neither predicted by a diffu- 
sion model, nor by a pure convection model. In our 
own Galaxy, this effect is more difficult to be clearly 
established |105 , and its effects have been investi- 



gated in particular on nuclei in various models |106| 



Tab. 2.4.9 compiles values obtained in several stud- 
ies for very simple forms of galactic wind. As one 
can see, almost all the authors set only upper limits 
on this possible wind. Actually, the exact form of 
galactic winds is not known. From a self-consistent 
analytical description including MHD calculations 
of the galactic wind flow, cosmic ray pressure and 
the thermal gas in a rotating Galaxy, a wind in- 
creasing linearly with z up to z ^ 15 kpc, with a 
z = value of about 22.5 km s~^ was found in |^ 
(see also references t here in). Following a completely 
different approach, |107| extract the velocity form 
able to reproduce data from a one-dimensional diffu- 
sion/convection model. They obtain a decrease from 
35 km s^^ to 12 km s^^ for z ranging from 40 pc 
to 1 kpc followed by an increase to 20 km s~^ at 
about 3 kpc. These two recent works raised a pos- 
itive conclusion about presence of a wind. For ref- 
erence, the values for our best B/C fits correspond 
to about 15 km s~^. The difficulty to compare con- 
stant wind values to other z-dependences is related 
to the fact that cosmic rays do not spend the same 
amount of time at all z, so that there cannot be a 
simple correspondence from one model to another. 
As a result, all the above-mentioned models are for- 
mally different, with, moreover, many different in- 
puts (spectral index, diffusion slope) that complicate 
a possible comparison. Nevertheless, their values are 
roughly compatible, the grounds for a physical mo- 
tivation for this wind being provided in . One can 
notice that other models (e.g. (2^) dislike convec- 
tive wind in their fit, but we will come back to this 



point in Sec. 5.3, To conclude about galactic winds, 
a radial dependence is not excluded Eq] . 



Reacceleration The term acceleration is used for 
the process by which low energy nuclei gain a huge 
amount of energy in a short time, promoting them 
to the rank of Cosmic Ray Nuclei. Once the cos- 
mic rays are accelerated, there are other processes 
that can lead to less sudden energy gains, referred 
to as reacceleration. There arc two classes of such 
proce sses. For the first, called sporadic reaccelera- 
tion 1 108 1, it is assumed that there are some well 
localized reaccelerating centers (e.g. supernova), on 
which the cosmic rays can scatter and gain a small 
amount of energy at each of these sc atter ings. For 
the second, referred to as continuous |10£|, reaccel- 
eration occurs during the wandering of the charged 



quence of spatial diffusion. The quasi-linear theory 
of Boltzmann equation, where all effects are aver- 
aged oyer statistical properties of the plasma (see 
e.g. [110|), indicates that the net effect is a diffusion 



in energy space 1 1 1 1 ] . The typical times for diffusion 
in space and energy are related through the relation 
Tspatiai-Trcac « L"^ /V^ (scc e.g. [|ll|l) , whcrc K IS the 
Alfvcn speed of the medium. Some basic features 
of these two modes of reacceleration can be found 
in Giler et al. (1989) pOSl. It is likely that both 



these modes are present (see as a good example, the 
treatment of |113|). 



The most direct observational evidence for reac- 
celeration comes from cosmic ray species which are 
unstable via electronic capture. Their abundance rel- 
ative to stable isotopes shows that they had some 
time to decay and therefore to attach an electron, 
which is more difficult at high energy. This means 
that these EC unstable species had a lower energy at 
the beginning of their propagation stage: the nuclei 
have gained about 100-200 MeV/nuc for a kinetic 
energy of a few hundreds of MeV/nuc [114|. A less 
direct evidence is related to the question of isotropy 
of cosmic rays measurements, which favors small dif- 
fusion slope (e.g. Kolmogorov S = 1/3) whereas the 
classical propagation models (no reacceleration) pre- 
fer larger values S ~ 0.6 (see and [ [109| ). It turns 
out that reacceleration allows to redeem small 6. 

In our studies |2^, ^ , we further investigated the 
allowed values of the reacceleration parameter, i.e. 
the Alfven speed Va. We tried to compare these val- 
ues with those used by other authors ^ ^ . The 
difhculty is that the secondary to primary ratios are 
determined by an effective value: 



ytruc ^ 



(18) 



where h and /iroac are respectively the height of the 
diffusive and reacceleration zone. The parameter w, 
which may depend on z, characterizes the level of 
turbulence and is often set to 1 [|l|. Our model, as 
others, uses 



uj{z) 



1 if Z < /iroac 

otherwise ; 



(19) 



as a crude approximation of the more complex real- 
ity. The total reacceleration rate, at least in a first 
approximation, is given by the convolution of the 
time spent in the reacceleration zone with the corre- 
sponding true Alfven speed in this zone. There are 
no direct observational clues about the size of the 
reacceleration zone, or about (jj{z). As a result, even 
if the analysis of cosmic rays can give values of ^ 
this does only give the value V2'"° up to an unknown 
factor h/hrcac- The situation is even more complex if 
no assumption is made about w(z). If w(z) strongly 
depends on z in a large reacceleration zone, then the 
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^-dependence 



model 



resolution 



data 



ref. 



^ 60 
- 8 
^ 20 

< 16 
^ 20 

< 15 
^ 20 
^ 20 

< 7 



constant (Vc) 
linear {2Vcz) 

constant (Vc) 

linear {3Vcz) 
constant (Vc) 

linear {Vcz) 



thin disc+halo (ID) 



thin disc +halo (2D) 



Monte-Carlo 
analytical 
numerical 
analytical 



numerical 



B/C, i"Be/"Be 



i"Bc 
B/C, i"Be/^Be 



(Owens, 1977) 
(Jones, 1979) 
(Kota & Owens, 1980) 
(Freedman et ai, 1980) 

(0, 1992) 
(Bloemen et ai, 1993) 
(Lukasiak et ai, 1994) 
(ItI, 1998) 
([271, 1998) 



Table 5: Constraints on galactic convective wind obtained from nuclei adjustments in various models listed 
in 



an effective size depending on the vertical distribu- 
tion of cosmic rays. However, this situation seems to 
be unfavored by MHD simulations . 

To conclude with an order of magnitude, our best 
models indicate that h/h^cac must be '--^ 1/4 in or- 
der to give realistic values |^ for Va (with lo ^ 1) p^ . 



2.5 Different approaches of propagation 

The equations describing diffusion have been writ- 
ten down in Sec. 2.3.(:, as Eq. (^). There are many 
approaches to solve the same diffusion equation, i.e. 
direct resolution (if possible) , extraction of the Green 
function of the problem, or separation of astrophys- 
ical and nuclear aspects through the weighted slab 
technique. Most of the theoretical work on cosmic 
rays makes an extensive use of the weighted slab 
model and of the so-called Leaky Box model. We 
now present a synthetic overview of these models and 
give some clues about their successes and failures. 



2.5.1 Leaky Box 

One of the simplest (though very simple) model is 
the so-called Leaky Box, in which the Galaxy is de- 
scribed as a finite propagation volume, delimited by 
a surface. Inside this volume, the densities of sources, 
interstellar matter and cosmic rays are homogeneous. 
Moreover, each nucleus has a probability per unit 
time 1/Tesc to escape from the box. In the station- 
ary regime, the densities are given by 

— + nva^N'^ = g^' + ^ nva^^ N'' . (20) 

'^esc , . , 

heavier k 

This model has been successful to explain most of 
observed cosmic ray stable fluxes at different ener- 
gies by a single function Tesc{E). This function can 
be either adjusted to the data, its physical interpre- 
tation being found afterwards, or extracted directly 
from more complete propagation equations [ST) . The 

3 There is a misnrint in where 1/4 .should be read in- 



second approach is a good way to understand why 
Leaky Boxes work so well, and is more easily seen 
in the framework of the weighted slab formulation. 
Thus we now turn to this point, introducing first the 
slab model. 



2.5.2 Slab model 

Independently of the exact processes which are re- 
sponsible for propagation, the first certainty we have 
about cosmic ray nuclei is that they have crossed 
some interstellar matter between their creation in the 
sources and their detection. This leads to spallation 
reactions that alter the initial (source) composition 
by destroying primary species and producing secon- 
daries. The ratio of secondary to primary species 
fluxes gives information about the quantity of mat- 
ter crossed. It is convenient to introduce the column 
density of matter crossed by a particle, also called 
the grammage, expressed in g cm^^. All the nu- 
clei of a given species, with a given energy, do not 
have the same propagation history, in particular they 
have not crossed the same amount of matter, so that 
a distribution of grammages is associated with each 
species. 

In a first step, though, one can assume that all nu- 
clei of a given species have crossed the same gram- 
mage. This is called the slab model. Formally, the 
number (x) of nuclei j that have crossed the gram- 
mage X is related to the destruction rate (the inelas- 
tic cross section is denoted ) and the creation rate 
(i.e. the spallation rate of all heavier nuclei k giving 
j, the cross section being denoted a'^^) by 



ax m 



-N'^ix) , (21) 



heavier k 



with the initial condition A^-'(a; = 0) = . The 
resolution of this equation yields the secondary to 
primary ratios as functions of x. Comparison with 
observations give the value of x. This leads to a con- 
tradiction, as weighting correctly the production rate 
by all the parent nuclei, the observed LiBeB/CNO 
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sub-Fe/Fe~ 1.5 gives x ^ 0.8 g cm ^ (see e.g. [ 115 



2.5.3 Weighted slab and other propagation models 

The previous model is too simple and one should take 
into account that the nuclei of the same species and 
same energy may cross different amounts of matter. 
Introducing the probability G{x) that a nucleus j has 
crossed the grammagc x, the density of a nucleus j 
is given by 



N^x)G{x)dx 



(22) 



The function G{x) is called the Path Length Distri- 
bution (hereafter referred to as PLD). It is remark- 
able that the use of the same function G{x) for all 
nuclei can solve the apparent contradiction between 
the grammagc seen by CNO and Fe: the Fe nuclei 
have a larger destruction cross section, so that they 
are more easily destroyed. The associated function 
N^{x) decreases more rapidly and the sub-Fe/Fe ra- 
tio is more sensitive to the low x part of the gram- 
mage distribution function G{x). 



ing feature of this model is that the average gram- 
mage, given by (x) = J xG[x)dx is exactly the es- 
cape length, i.e (x) = Xcsc{E). Actually, this quan- 
tity solely depends on energy if G^^{x) is the same 
for all species. This PLD may be modified to give 
new weighted slab models. For example, some dis- 
crepancies have led authors in the past to introduce 
the double PLD |118| which is a mixture of two 
G^^{x) with two different characteristic functions 
Aosc and X'^^^^. This can be interpreted as the pres- 
ence of two different propagation zones (Simon et at, 
1977, 1979 |11S|), which supports the more realistic 
disc-halo model. Let us notice that the terms escape 
length and Path Length Distribution are sometimes 
mixed in the literature, which can generate confu- 
sion. 

Finally, it is worth noting a specific point quite in- 
dependent of th e rest of the discussion. It has been 
noticed in 119 that a different but equivalent for- 
mulation could be obtained from the above set (23) 
of LB Weighted Slab representation, by the substi- 
tution 



GLB(a;) = (l/Aesc)cxp(-x/Aesc) ^ G{x) = 1 



Weighted slab as models or as a general tech- 
nique The probability distribution G{x) can be de- 
termined by the choice of a propagation model. For 
example, the slab model is obtained by considering 
G{x) = S{x — xq). Actually, the weighted slab ap- 
proach has been widely used in literature, under two 
slightly different forms. In the weighted slab model, 
one tries empirically to modify the Path Length Dis- 
tribution G{x) to account for the data. The phys- 
ical meaning of G{x) may be explored in a second 
step. The weighted slab technique is more general, 
and the name refers to the possibility to introduce 
Eq. ( p2[ ) containing the additional and meaningful 
variable x (the grammagc) in any propagation equa- 
tion (such as the LB model, diffusion model, diffu- 
sion/convection model...). Exa mple of d irect ex- 
traction of PLD can be found in [116, |ll7|[ . Up to 



the end of this section, we will omit this subtility. 
The weighted slab approach allows to link Leaky Box 
models with more realistic diffusion models, explain- 
ing why these Leaky Boxes work so well. 

Path Length Distribution of the Leaky Box 
model As an alternative to the direct analytical 
resolution of Eq. (^), one can insert Eq. (|2^ ) in 
Eq. (|20| ) which separates the nuclear part from the 
other effects. This gives 



G^^ix) = ^exp(-^ 

^csc \ -^csc 



dx 



a-' ~ ■ — ^ cr^-^ ~ 1 

m ' ^ m 



(23) 



and the addition of a term N^/X^ 
part of Eq. (H), i.e. 



in the nuclear 



dNHx) 
dx 



kj 



- — + —N^ix)=y —N^ix). 

k>3 



The meaning of the new set of equations is clear: ei- 
ther one solves the stationary Eq. ( pO| ) , or one solves 
the same but non- stationary equation integrating 
over aU times, i.e. / N^{x)G{x)dx = J N^{x)dx 
(here, cc is a dummy variable, so that it can be called 
time). Such a procedure is very general. It is for 
example used in [120[. 



Weighted slab technique applied to diffusion 
models The link between the diffusion equation 
and the escape time introduced in Leaky Box mod- 
els was clarified by Jones in [121 . He showed that 



k>j 



diffusion models, energy losses included, can be rein- 
terpreted in terms of a Leaky Box model (with some 
exceptions, such as e~ for which the synchrotron or 
inverse Compton losses arc too important). This 
reinterpretation, called the leakage lifetime approx- 
imation can be understood by writing the diffusion 
equation, neglecting energy gains and losses and as- 
suming steady state, 

-V.{K^VN^)+nv<7^N^ = q^^nva'^^N'' .(24) 

k>j 



In this expression, the PLD G(x) is expressed as 



When applying the weighted slab technique (Eq. p^) , 
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fusive and a purely spallative part 



-/vAG(r,a;) = Q{r,x) ; 

' —N^{x)^} —N^{x) 
m ^-^ m 

k>j 



dx 



(25) 



The Path Length Distribution G(r, x) encodes all the 
propagation properties. This function has been ex- 
tracted for various geometries (thin disc, spherical 
halo) and various forms of source spatial distribution 
Q{r) |116|. A general result is that as the operator A 
is hermitic, the function G can always be set under 
the form 



G(r, x) (X ^ c„ (r) exp f - ^ 



Comparing to G^^ {x) derived above in Eq. ( p3| ) , this 
shows that diffusion models can be equivalently ex- 
pressed in terms of a sum of Leaky Boxes. Jones 
showed 1 121 that in most situations, only the first 
terms contribute in the sum above. This explains 
why the simplest Leaky Box models are so success- 
ful to describe diffusion, even in a complex geometry. 

The specific case of primaries is noticeable. In this 
case, the nuclear part of Eq. ( ^ ) reduces to a very 
simple form to finally give 



N\r) ~ / G(r,x)exp{—a^x)dx 



(26) 



Mathematically, the solution is the Laplace trans- 
form of parameter cr^ of the function G{r,x). Sev- 
eral solutions corresponding to several diffusive ge- 
ometries can be tabulated with this Laplace trans- 
form 1 122]. This approach gives an alternative way 
to evaluate the average grammage (a;), i.e. the equiv- 
alent Leaky Box description, from 



X G{x)dx 
G{x)dx 



da^ 



To end with these relations between treatment of 
the same equations, the latter formulation is close 
to what one obtains using the random walk ap- 
proac h ||65| , 123 1. This should be not too surprising 
since [|l7| showed a long time ago the equivalence of 
random walk description and diffusion processes. 



2.5.4 Limitation and usefulness of all these models 

Actually, the fluxes derived with the weighted slab 
technique as depicted above differ slightly from the 
direct evaluation. Whereas Lezniak 124 tried to 
extend the simple scheme p rese nted in Eq. (|2^) to 
include energy losses, Jones [125| demonstrated that 
the separation between the nuclear side and the prop- 



high energy regime. The reason is that in the rewrit- 
ing N^{r) = N^{x)G{r,x)dx, the separation be- 
tween (x) and G(r, x) can have only one param- 
eter in common, otherwise, equations obtained from 
the initial propagation equation are generally incon- 
sistent. Hence, one can think that energy appears 
only as a mute parameter, but this is not exactly the 
case: the astrophysical part, i.e. G(r,x) is the same 
for all species as long as rigidity is considered. Con- 
versely, for the nuclear part, the equation obtained 
above is only valid for a fixed kinetic energy per nu- 
cleus. 

However, redefining some parameters and under 
seve ral c onditions, the technique can be made ex- 
act [126 , and this idea ha s be en numerically and 
quantitatively validated in |12C|. Despite that, the 



leakage lifetime approximation is not valid in some 
cases a mentioned above. For some nuclei, a descrip- 
tion in terms of Leaky Box may lead to wrong results. 
This is the case for radioactive species in a realistic 
Galaxy, as was shown by |127|. 



Concluding remarks The Leaky Box models, 
due to their simplicity, arc very well suited to the 
extraction of source abundances (elemental as well 
as isotopic). They can also be used to compute 
the secondary antiproton production, since the same 
processes as for secondary stable nuclei are at work. 
However, as emphasized in (see also Sec. 2.4.6), 
they are not able to predict any primary contribu- 
tion in the antiproton signal, since it requires the 
knowledge of the spatial distribution of primary pro- 
genitors. It is also well known that the Leaky Box 
parameters are just phcnomcnological with only a 
distant connection to physical quantities. However, 
even if apparently many existin g eff ects cannot be 
correctly included, it is shown in |l28[ that the galac- 
tic wind can be accounted for in a Leaky Box de- 
scription and emphasized that the phenomenological 
behavior of the escape length at low energy could be 
due to the presence of this galactic wind. This idea 
was investigated further in ||57|| with a generation of 
several equivalent phenomenological escape lengths 
from several possible physical configurations of a one- 
dimensional diffusion model. The relation between 
one-dimensional models and Leaky Box models is 
thus firmly established and very well understood. 
This relation also elucidates some of the physical con- 
tents of Leaky Box models. 

Our model furnished the following step towards 
realistic description, because it is equivalent, up to 
several minor modifications, to the one-dimensional 
model of (st) . Finally, in the Strong et al. model , 
all subtle effects can be studied and modelled, with 
the counterpart that the numerical approach makes 
the physical intuition of the results less straightfor- 
ward and the computation slower. However, the 
physical input is almost the same as in our model, 
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Even if all these models are equivalent to describe 
the local observations of charged cosmic rays, they 
lead to very different conclusions and interpretations 
when the spatial variation of the cosmic ray density 
is considered. As an illustration of the poor current 
understanding of this global aspect, we mention the 
ever-lasting problem of the gamma ray excess about 



1 GeV towards the galactic center |129| or the too flat 
radial 7-ray distribution observed in the disc [130|. 



2.6 Numerical implementation 

From the above discussions, it should be clear that 
the evaluation of the fluxes involves the computa- 
tion of a nuclear reactions grid and the resolution of 
differential equations describing diffusion. 



Nuclear part The flux of a given nucleus is the 
sum of the primary (source) contribution and all 
the secondary contributions, from the spallations of 
heavier nuclei. It may be convenient to write a 
ma trix relation (easily diagonalized, see coefficients 
in (m)) between the fluxes {Ari}J=i - " ^^^^ ^^le 
sources - " 

[N] = [aM + mN] , 

where the matrices [a] and [/?] contain the informa- 
tion about destructive spallations and secondary cre- 
ation. This is particularly useful in the Leaky Box 
model, as this relation is straightforwardly inverted 
and can provide source abundances as we ll as associ- 
ated errors without too much efforts 122 1 . How- 
ever, as soon as energy losses are considered, this 
becomes more complicated. 



We chose a more direct method, i.e. the cascade 
method. The nuclei are classified according to their 
masses. The heaviest nucleus A^i has no secondary 
contribution and is given directly by its source contri- 
bution qi. Eq. (H) is solved for this first nucleus. The 
flux of the second heavier nucleus N2 is a mere com- 
bination of the secondary contribution coming from 
A^i and its own source contribution 52- Then, Eq. (||) 
is solved for this second nucleus, using the result ob- 
tained for A^i . The flux of the third N3 is a combina- 
tion of its own source contribution (73 , plus spallative 
contributions from the two previous ones A'^i and 7^2- 
This procedure is repeated for all the nuclei. These 
implicit contributions are displayed in Fig. ^: source 
abundances have been arbitrarily normalized to 100 



(dashed right bars, Solar abundances |132| times first 
ionization potential, see e.g. |133 ). 



Left bars are for propagated fluxes at about 1 
GeV/nuc in a simple Leaky Box model without en- 
ergy losses. A first point is to notice that heavy 
nuclei are more subject to destruction during propa- 
gation than light nuclei. Secondary species also ap- 
pear clearly on the plot. Coming back to the cascade 
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[KKl Source abundances (before propagation) 

All secondary contributions added (i'uU sbower) 
I I Direct secondary contribution from primaries added 
I I Propagated source abundances (primaries) 
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Nuclei Z = 4-28 



Figure 9: This figure shows the source abundances, 
along with propagated abundances. For the sec- 
ondary species, we indicate the abundance that 
would be obtained by neglecting the secondary con- 
tribution due to the spallation of heavier secondaries. 
This second order effect is discussed in the text. For 
the sake of clarity, the scale has been expanded for 
species heavier than P, by a factor x20 for species 
lighter than Mn and x5 for heavier species. 



(left bars) the sole propagation of sources (dashed 
bars), i.e. primaries; faint shaded bars show the 
sole secondary contribution (direct contribution of 
heavier sources); the higher bars correspond to the 
summation of all intermediary steps. The effect of 
these higher order secondaries is particularly impor- 
tant for heavier species. An iterative method can also 
be used pTf , implementing the previous step propa- 
gated abundances as sources for the next step, until 
convergence is reached. 



Spatial part The spatial part of the diffusion 
equation involves second order differential equations. 
In our semi-analytical approach, one deals with a sec- 
ond order equation in energy but no spatial deriva- 
tive. The resolution is much simplified, compared 
for example to where the full second order par- 
tial derivative transport equation is solved. As re- 
gards the energy part, the Rimge-Kutta method 
is not suited and one must turn to more refined 
methods, such as the Crank-Nicholson scheme. The 
general method can be found in many numerical 
books 1 134 1 . An alternative approach is random- 
walk Monte Carlo simulations. A diffusive pro- 
cess can be mimicked by a random walk, with step 
related to the diffusion coefficient 0. This ap- 
proach is very time consuming but a few authors has 
used th is te chnique since mid-seventies (e.g. Owens, 
1976b 1 124]) to p ropagate self-consistently proton 
and electrons |l35|] . Such an approach is useful when 
one wants to model very complex and inhomoge- 
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2.7 Quick survey of charged nuclei behaviors 
and their interest 

We now have all the elements needed to compute 
the fluxes of all the cosmic ray nuclei. They all have 
different propagation histories, and they carry infor- 
mation of different nature (sec [136| for a review). 
They can be classified as follows: 

• Stable primaries and mixed species (i.e. 

secondary plus primary contributions): most 
elemental abundances point towards a "stan- 
dard" origin (from a particle physicist point of 
view). Some nuclei present anomalous contri- 
butions that could be due to wind-induced en- 
richment from Wolf-Rayet stars surfaces. The 
correct abundance determination of nuclei be- 
longing to this class also allows a better charac- 
terization of the acceleration processes. 



• Pure stable secondary species: the three 
groups of secondary species {Z = 1 — 2, Z = 3—5 
and Z = 21 — 23) give information on the propa- 
gation history for quite different charges. Their 
importance for the determination of the prop- 
agation parameters is extensively demonstrated 
in this review. 



P unstable secondary species: though ra- 
dioactive species were first used to show the ex- 
istence of a large diffusive halo [137|, we under- 
lined in [|0| the interest of ^°Be, ^^A\ and ^^Cl 
(^^Mn is more difficult to tackle since it has a 
mixed decay mode, f3 and EC) for the study the 
LISM, or as a possibility to derive the diffusion 
coefficient once LISM properties are fixed [pl[. 



• Primary EC clocks: through Co/Ni, Co/Fe 
or the more precise ^^Ni/^'^Ni ratios, these EC 
clocks indicate the time elapsed between syn- 
thesis and acceleration and are shown as circles 
in Fig. ^ (the latest ACE results point towards 
< > 5 X 10^5 yr 111). 



• Secondary reacceleration indicators: other 
K-capture nuclei (open diamonds in Fig 



loD 



indicate that reacceleration during propagation 
could be as large as a few hundreds of MeV/nuc 
at low energy. The latest results from the ACE 
experiment seem to confirm this reacceleration, 
but the propagation of EC unstable nuclei in a 
refined propagation model and a thorough anal- 
ysis including peculiarities of the LISM has not 
yet been performed (Maurin et a/., in prepara- 
tion). 
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Figure 10: Stable and unstable cosmic rays in prop- 
agation model and their primary or secondary char- 
acter (taken from |138|). 



3 Solutions of the diffu- 
sion/convection equation for two- 
zone cylindrical models 

3.1 General remarks 

The density of cosmic rays of energy E at the 
position (r, z) is obtained by solving the energy- 
dependent diffusion equation (see Sec. 2.3.6 ). It 
turns out that it is possible to first focus on the spa- 
tial diffusion and, then to take the energy changes 
into account, as long as these energy changes take 
place in the disc. The density is then obtained by 
solving a Laplace equation in a cylindrical geometry. 
The standard method is to develop all the quantities 
over a suitable complete set of orthogonal functions 
involving the first Bessel function Jq. Introducing 
p — r/ R, one can write 

oo 

N{r,z) = Y,N,{z)UQp) , 



i=l 



with 



N^{Z) 



JiiQ) Jo 



pN{pR,z)Jo{Cip)dp 



the Q are the successive zeros of the function Jg. 
This development is inserted into the diffusion equa- 
tion, which can then be rewritten, using the proper- 
ties of Bessel functions, as 



Ni'iz)~(^^Siz) + ^ + ^]miz) 



■^N^iz) 



= (z) 



K 



(27) 



where Finei and Frad denote respectively the spal- 
lation rate over the interstellar medium and the ra- 
dioactive decay rate. The Bessel expansion qf°^'''^°{z) 
for a unspecified source term q(r, z) has also been 
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3.2 Full solutions 

We do not recall here the procedure used to solve 
Eq. (H) and we refer the suspicious reader to p9| . 
Suffice to say that very general solutions can be writ- 
ten for each species. They consist in the sum of two 
terms, the first involving sources located in the disc 
(primary or secondary), and the second involving 
sources located in the whole halo (radioactive ori- 
gin). For the sake of clarity, we present separately 
these two contributions. 



3.2.1 Contribution from progenitors in the disc 

The galactic disc is a particular place for GCR, as 
it contains the standard primary sources q{r, z) = 
(7'*""^(r)5(z) (sec Sec. 2.4.3 ), and it is the place where 
spallations occur. The contribution to the density is 
given by 



ATI ^ (^-^ 

N{r,z) = cxp ( — 



E 



smh[S,iL-z)/2] 



^ A, smh[S,L/2] 



(^) (28) 



with 



2hr 



incl 



Vc + KS, coth 



i?2 ' A'2 ^ K 



S,.L 



For a pure primary Qj = qQQ{E) x qf^^'^ - where 
^disc denotes the Bessel expansion of q{r) - and for a 
pure secondary Qi = ^™''>™3 T''^N,^{0) ~ where the 
superscript j denotes the nucleus evaluated, omitted 
in above expressions -. 

When energy losses and diffusive reacceleration are 
taken into account, the differential equation is ob- 
tained, following the procedure described e.g. in p9| , 
as 



^"Be— ^^''B, neglecting this channel would give an er- 
ror of about 10% on the B flux, whereas considering 
that this term is only located in the disc would give 
an error of about 3% compared to the rigorous treat- 
ment given above. Due to lack of space, we do not 
reproduce the expression in this case, but the reader 
is referred to [|9| . 

3.2.3 Contribution from the halo: primary sources 

We will consider the hypothesis that primary sources 
of cosmic rays may be present in the diffusive halo 
(see Sec. ^.4.4 ). The solution of the diffusion equa- 
tion is somewhat different than in the standard case 
of sources located in the disc. Denoting ^'^^1° the 
corresponding source term, the Bessel terms of the 
density are given by (see |139[|) 



N^iz) 



y,{z) y,(L)e^-(i-l-^)/2^ 
' KS, A,sinh(S',;L/2) 

Vc + 2/iri„el 



cosh(S'iz/2) 



sinh(S'iz/2) 



with 



yi{z) = 2 1 cxp 
/o 



Vc{z - z') 



X sinh 



2K 

Si{z — z')\ halo/_/\ 



M halo/ /\ 7 / 

-jq, {z )dz 



In particular, the density in the disc {z = 0) is given 

by 



7V,(0) = exp 



VcL 



2K J Aifi\nh{S,L/2) 



It can be checked that with a source term localized 
in the disc, i.e. g*^'^'°(r,z) = 2M(z)g'i'"'=(r), the ex- 
pression (^) is recovered. 

3.2.4 Modelling of the local bubble 



d 



(29) 



dE 



d 



2h— btot{E)N,{0) - p'Kpp—N,{0) 



dE 



where the energy loss a nd r eacce leration terms were 
discussed in Sec. 2.3.3 and 2.3.5. Subscript i refers 



to Fourier-Bessel coefficients for the above equation 
expended on cylindrical Bessel basis Jo(C,ip). 

3.2.2 Contribution from the halo: radioactive /3 decay 

All the nuclei treated here have at most one unsta- 
ble progenitor. The contribution of these radioac- 
tive nuclei may be unimportant in some cases, but 
we should take it into account as it is the dominant 



As mentioned above, it is very likely that the Sun 
lies in an undcrdense bubble. The main effect on 
the cosmic rays is that spallations, and hence the 
secondary production and destruction, are locally 
strongly reduced. In the thin disc model, this can 
be taken into account assuming position dependent 
spallation terms. Actually, as will be discussed below 



(see Sec. 6.3), this effect is not so important for stable 
species, but it is crucial for the low energy radioactive 
species, which are only sensitive to the local propaga- 
tion conditions. Once created at a given point, they 
diffuse in a small region surrounding this point before 
they decay. The density of their progenitors does not 
vary much in this region, and it can be considered 
as homogeneous. Moreover, the propagation of the 
radioactive is almost not altered by the boundaries 
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center of the diffusive volume at the position of the 
Sun. The bubble can then be modelled as a hole in 
the disc. The rates TspaK^) nusuir) .iKaspa^i and 
rinei(r-) = ni^isMir).v.<7inci now depend explicitly on 
r via the local interstellar density which reads 



riLiSMir) = e(r - Hioio) "ISM 



(30) 



where Q is the Heaviside distribution (sec Fig. pi] ). 
As for the no-hole model, a solution is found by 




Figure 11: Schematic representation of the hole 
model. The disc has zero thickness with a hole of 
size Thole in gas density. 

expanding the density over the Bessel functions 
Jo{(,iT/R)- The expression for the solution may be 
found in iol. 



3.3 Physical interpretation of the formulae 

For the sake of clarity, we wish to further discuss two 
particular cases. To begin with, the case of primary 
stable species is of particular interest. Their density 
in the disc is given by 



„disc 



which, when compared to the Bessel expansion of the 



q{r) 



contains basically all the physics of diffusion. Diffu- 
sion in Bessel space may be seen as a filtering process, 
with a filter 1/Ai. The small scale features, which 
correspond to large values of the index i, are spread 
out and erased by diffusion as Ai 2K(i ^ irKi for 
large i. The large scale structure is determined by 
the small i behavior. The effect of small L, spalla- 
tions and galactic wind is to limit the decrease of Ai 
for small values of i, i.e. they filter the large scale 
features out. In other words, they make diffusion a 
local process. More precisely for small L and small 
i we have Ai « 2hTiac\ + Vc + 2K/L. It appears 
that three characteristic lengths, although of differ- 
ent origin, play a very similar role in erasing the large 
scale features of the source distribution, namely L, 
Twind = 2K/Vc and rgpai = 2K/2hTiaci- Structures 
larger than these are destroyed by diffusion. This 



4 Origin of cosmic rays 

4.1 The question of the spatial origin 

The cosmic rays emitted from distant sources are 
more likely to escape from the diffusive halo or to be 
destroyed by spallation before reaching Earth than 
those coming from nearby sources. Conversely, a CR 
eventually detected on Earth has a greater proba- 
bility to come from a nearby source than a distant 
source. In this section, we investigate the question 
of the spatial origin in more details, and we find that 
under quite reasonable conditions, most of the cos- 
mic rays detected in Earth were emitted by sources 
located in a radius less than ~ 3L. 



4.2 Method 

The question we wish to address is the following: a 
cosmic ray being detected at the position Tq of an 
observer (in practice, this will be the position of the 
Sun), what is the probability 



dV {emitted : rs,^s + dfs 
= V{rs\ro} drs 



observed : Vo} 



that it was emitted from a source located at the po- 
sition Ts? Such a question falls among classical prob- 
lems of statistics. A rigorous theoretical frame is pro- 
vided by the Bayes approach that summarizes proper 
use of conditional probabilities. A cruder but suffi- 
cient (and equivalent) treatment is given through the 
frequency interpretation. The probability is simply 
given by the number of particular realizations, i.e. 
paths M [ts To] (the cosmic rays travels from the 
source at position Ts to the observer at position fo) 
divided by the number of all realizations of Af [-^ fo] 
(the cosmic rays reach the observer position). No- 
tice that the latter number is just a normalization 
so that the first quantity determines completely the 
problem. In the frequency interpretation, we thus 
have 



r{fs\fo} = 



(31) 



We finally notice that the number of paths 
Af [rs — ^ To] determines the number of cosmic rays 
that reach position fo, when a source is placed at 
position fs- We can thus write 



r {fs\fo} 'X N.^ifo) 



(32) 



where the density N,:^{fo) is the solution of the prop- 
agation equation for a point source located at r*,. 
The normalization factor of this relation is obtained 
by imposing that V actually is a probability, i.e. is 
normalized to unity. If the sources are distributed ac- 
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detected at was emitted from a volume V' (or a 
disc surface) is given by 



V{V\ro} = 



(33) 



This probability contains all the physical information 
about the spatial origin of cosmic rays and may be 
computed for different situations (sec for further 
details). 



4.3 Pure diffusion case and key parameters for 
spallations and wind 

The extent of the zone from which the cosmic rays 
detected on Earth actually come from is limited by 
all the phenomena that keep them from propagating 
over large distances: escape through the boundaries, 
spallations and galactic wind. It comes as no sur- 
prise, then, that we can define three typical lengths, 
beyond which the cosmic rays do not reach us 



rL 



T'spal = K/hTinc\ 

w 3.17 kpc X 

''wind = 2A7K 

w 5.87 kpc X 



(34) 



100 mb 



0.03 kpc^ Myr^^ 



K 



0.03 kpc^ Myr"^ 



(35) 

10 km s-i 



The role of the length r^ind hi the context of spatial 



origin was recognized in |128|. For a given set of dif- 
fusion parameters, these lengths can be computed, 
as well as the map of the probability that a particle 
reaching Earth was emitted from each source located 
in the disc. The lengths r-wind and rgpai depend on 
K , and are larger at higher energy. This is of course 
not the case for = L, so that at sufficiently high 
energy, r^ind > L and rgpai > L, the spatial origin 
being solely dictated by the halo size L (side bound- 
ary r = i? is unimportant). The purely diffusive case 
is illustrated in Fig. |l^, which shows the regions from 
which a given fraction of the cosmic rays reaching the 
Earth were emitted from. We refer the reader to 
for the details. 

Conversely, it turns out that when realistic diffu- 
sion parameters are considered, these lengths may be 
rather small, in particular for heavy nuclei which are 
more sensitive to spallations, so that the observation 
of cosmic rays in the Solar neighborhood only gives 
information about the local conditions under which 
diffusion occurs. This will be discussed in more de- 
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Figure 12: This figure shows the regions of the disc 
from which a given fraction (50 % for the inner sur- 
face and 99 % for the outer surface) of the cosmic rays 
reaching the Earth were emitted from. The dashed 
lines represent the case of a halo with no upper and 
lower boundaries {L oo) whereas the solid line 
are for L = 2 kpc. The spallations and galactic wind 
are not taken into account here, so that these regions 
overestimate the origin of the sources (they represent 
the geometrical limit). 

5 Experimental determination of the 
diffusion parameters 

We wish to study the constraints on the diffusion 
parameters (the spectral index of sources a, the nor- 
malization Kq and the spectral index d of the diffu- 
sion coefficient, the height of the diffusive halo L, the 
galactic convective wind speed Vc and the Alfvenic 
speed Va) that come from the measured fluxes of cos- 
mic rays. For the aim of our analysis, we can consider 
different classes of flux ratios: primary-to-primary 
(e.g. C/0), secondary-to-primary (e.g. B/C or sub- 
Fe/Fe), secondary-to-secondary (e.g. Li/B or Be/B), 
ratios of either stable (e.g. ^"B/^^B) or unstable (e.g. 
-'^"Be/^Be) isotopes. Each of these may be an indica- 
tor of some dominant physical phenomenon and be 
particularly sensitive to the corresponding diffusion 
parameters (see also Sec. 2.7). The ratio of two pri- 
maries is practically insensitive to changes in all the 
parameters, since they have the same origin and un- 
dergo the same physical processes (but keep in mind 
that extreme nuclei such as p and Fe present dras- 
tically different destruction rates, see also Sec. 6.4). 
Similar conclusions, even if less strong, may be drawn 
for the ratio of two isotopes of the same species, such 
as ^"B/^^B. Indeed, at very low energy values this 
quantity is slightly affected by changes in the injec- 
tion spectra, but the effect is too weak to constrain 
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One of the most sensitive quantity is B/C, as B 
is purely secondary and its main progenitors C and 
O are primaries. The shape of this ratio is seriously 
modified by changes in the propagation coefficients. 
Moreover, it is also the quantity measured with the 
best accuracy, so that it is ideal to test models. In- 
deed, as a ratio of two nuclei with similar Z, it is 
less sensitive to systematic errors and to Solar mod- 
ulation than single fluxes or other ratios of nuclei 
with more distant charges. For the same reasons, the 
sub-Fe/Fe may also be useful. Unfortunately, since 
existing data are still affected by sizeable experimen- 
tal errors, we can only use them to cross-check the 
validity of B/C but not to further constrain the pa- 
rameters under scrutiny. 



5.1 The method 

For a given set of parameters, the source abundances 
of all nuclei (i.e. primaries and mixed nuclei) are 
adjusted so that the propagated top of atmosphere 
fluxes agree with the data at 10.6 GcV/nuc (sec Q). 
The nuclear cascade is started from Sulfur, as we 
checked that the heavier nuclei do not contribute sig- 
nificantly to the B/C ratio. The top of atmosphere 
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Figure 13: Diagrammatic representation of the vari- 
ous steps of the propagation code 

fluxes arc deduced from interstellar fluxes using the 
force field modulation scheme (see EgII and references 



pared to the data (see below) and the is computed 
for the chosen set of parameters, in our case 



E 



((B/C). 



(B/C)i, model) 



(36) 



2,exp 



where the sum runs over 26 experimental values from 
HEAO-3 [ |140[ . In general, if the experimental set-up 
is such that the measured (experimental) values dif- 
fer from the "real" values by a non biased quantity 
with a given probability distribution, then the 
value gives a quantitative estimate of the probability 
that the model is appropriate to describe the data. 
However, this condition is probably not fulfilled for 
HEAO-3, as for some measured quantity, the quoted 
errors crf^^^ are much smaller (e.g. oxygen fluxes) 
or much larger (e.g. sub-Fe/Fe ratio) than the dis- 
persion of data themselves. For this reason, it is 
meaningless to associate a likelihood to given val- 
ues. As a consequence, in an ideal situation in which 
very good and consistent data on B/C and sub-Fe/Fe 
ratios were available, the best attitude would be to 
make a statistical analysis of the combined set of 
data. Unfortunately, this is not currently the case. 
One can follow two ways to extract information from 
the sub-Fe/Fe data. First, as a check, one can com- 
pare the sub-Fe/Fe ratio predicted by our model - 
using the parameters derived from our above B/C 
analysis - with data from the same experiment. Sec- 
ond, one can search directly the minimum x^e 
sub-Fe/Fe ratio, with no prior coming from B/C. 
This procedure is more hazardous since the statis- 
tical significance of the sub-Fe/Fe data is far from 
clear. 



5.2 The available data sets 
5.2.1 B/C and sub-Fe/Fe 

As emphasized above, the results presented here are 
mostly based on the data taken by heao-3 |140| . 
They have been taken in 1979-80, for elements with 
charges from 4 to 28 and for energies ranging from 
0.6 to 35 GeV/nuc around a minimal Solar activ- 
ity. In the case of B/C, the quoted 1-cr HEAO-3 rel- 
ative errors are 2-3%. We also considered data from 



balloons |141] and from the isee-3 experiment ]142t , 
even if the relevant error bars are wider. The first one 
collected data in 1973-75 for energies spanning from 
around 1.7 to 7 GeV/nuc. The second experiment 
was operating during 1979-81 on board a spacecraft, 
in the energy range 100-200 MeV/nuc. In some of 
the figures presented below, we also plot - for purely 
illustrative goals - the data point from imp-8 [|| and 
the VOYAGER experiments |142] (we did not include 
ULYSSES data point |^ , since it corresponds to a pe- 
riod of maximal Solar activity). Nevertheless, we 
have checked that the best x^ values were not signif- 
icantly modified when these points were added. 
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used data from heao-3 140 and from balloons |141|, 
In both cases the error bars, around 10%, are signif- 
icantly larger than for B/C. 



5.2.2 Radioactive species 

Several experiments in the last twenty-thirty years 
have measured radioactive isotopes in cosmic rays 
with increasing precision, at energies of a few hun- 
dreds of MeV/nuc. The first data - usually presented 
as the ratio of some radioactive isotope to its sta- 
ble companion(s) - were affected by errors of around 
25-30%. The latest published data have error bars 
reduced by a factor of two or three. In the follow- 
ing, we implicitly refer to three satellite experiments, 
namely voyager, ulysses and ACE. Other experi- 
ments will sometimes be shown on figures but they 
will be purely illustrative since their accuracy is far 
smaller. 

The best measured ratio is probably ^"Be/^Be 
which corresponds to the lowest Z /3-radioactive nu- 
cleus. Data from ulysses |144| and from ACE [ |145| 
are consistent, the quoted error bars being smaller 
for ACE. They are also consistent with the voy- 
ager data point |143| for which the quoted error is 
larger. We do not use the SMILI data, as the possibil- 
ity that they are plagued by statistical fluctuations 
is not ruled out |146|. 

As regards the radioactive chlorine isotope '^^Cl, 
results are usually provided as '^^Cl to total CI ratio. 
The only available data, to our knowledge, are those 



from ULY SSES 1 147 with a l-cr error of about 35%, 
and ACE [145] whose errors (even taken at 3-(t) are 
completely included in the ULYSSES l-cr upper error 
band. 

Finally, the measurement of the ^^Al/^^Al ra- 
tio is m ore problemat ic. Indeed, the data from 
ULYSSES |14^ and ACE [145 1 do not seem to be com- 



patible (the ACE central point is much lower than 
ulysses' one). Even enlarging ACE error bars (which 
are smaller than ulysses) to 3-cr does not improve 
significantly the compatibility. On the other side, 
the ULYSSES data are fully compatible with 1-a VOY- 
AGER [149[ ones, whose uncertainty is still much 
greater than for the other two experiments. The 
possible discrepancy between some of these data is 
addressed in [50 . 



5.2.3 p and He fluxes 

We do not study here the compatibility of p and 
He fluxes with the prediction from diffusion mod- 
els. However, we need the observed values of these 
fluxes to compute the secondary antiproton fluxes. 
The contribution of heavier nuclei to the antiproton 
production is negligible. Until recently, the spectra 
of p and He were known with a modest accuracy 



incompatible at high energy. This induced an un- 
certainty of some tens of percents in the predicted 
antiproton spectrum. Recent measurements made 



by the balloon-borne spectrometer BESS [150[ and by 
the AMS detector during the space shuttle flight [151 
dramatically reduced the uncertainties both on pro- 
ton and helium spectra. We fitted the high energy 
(T > 20 GeV/nuc) part of these measured spectra 
with the power law: 



$(T) = N(T/ GcV/nuc)"'' 



(37) 



where the kinetic energy per nucleon T is given in 
units of GeV/nuc and the normalization factor N in 



units of m 



sr (GeV/nuc) . This provides 



a good description down to the threshold energy for 
the antiproton production. 

Wc fitted the bess and AMS data both separately 
and combined, obtaining very similar results. This 
is obvious since the data from the two experiments 
are now totally compatible. The best fit corresponds 
to iV = 13249 m-2 sr^^ (GeV/nuc)"^ and 

7 = 2.72; spectra obtained from the best fits on 
the single BESS and AMS data completely overlap. 
We did the same for helium and the corresponding 
numbers are N = 721 sr~^ (GeV/nuc)"""^ 

and 7 = 2.74. The 1-a deviation from the best fit 
spectrum docs not exceed 1% for both species. The 
uncertainty induced on the antiproton spectrum is 
smaller than the ones induced by the uncertainties 
on nuclear physics and diflfusion parameters. The 
situation has significantly improved since 1998, when 
an error of ± 25% was quoted [152[. 



5.3 The diffusion parameters estimated from 
HEAO-3 B/C measurements 

5.3.1 Partial exploration of the parameter space 

In a first step p9l , we did not consider a,- , the slope of 



the source spectrum of each species j (see Sec. 2.4.3 



p. |ll|), as a free parameter but instead we set it to 
the value obtained by subtracting S from the spec- 
tra measured at high energy (see [^). We display 
models that have < 40. 



Parameters are strongly degenerated Many 
configurations for L, Kq, S, Vc and Va lead to the 
same B/C ratio. This can be understood as meaning 
that they all lead to the same effective description 
in terms of the associated and grammage (see 
discussion in Sec. 2.5.3). In Fig. |l^, each value of d 



gives a different contour plot in the Kq/L — L plane. 
It appears that they all can be superimposed to a 
single curve by a rescaling Kq/L Kq/L x f{S), 
where / is a function of S only. For the contours dis- 
played in Fig. |lj, it takes the values /(0.46) = 0.51, 
/(0.5) = 0.62, /(Q.6) = 1, /(Q.7) = 1.54 and 



30 5 EXPERIMENTAL DETERMINATION OF THE DIFFUSION PARAMETERS 



the values of Vc are shifted downward as S is de- 
creased but the allowed range of Va/^/K^ does not 
significantly move. In particular, with the above cri- 
terion, we do not find any model having a good 
without convection {Vc = 0) or without rcaccelera- 
tion {Va = 0). We emphasize the fact that when a 
and 6 are fixed, the four remaining parameters are 
strongly correlated so that, when one of them is de- 
termined, in principle the allowed ranges for the oth- 
ers will be narrower than what could naively appear 
from the figures. 
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Figure 15: Models with different values of S, the dif- 
fusion coefficient spectral index, are shown. For each 
value of Vc and Vq/VKq , only the best value 
is retained when the other parameters L and Kq/L 
are varied. The figure displays the contour levels for 
< 40 for the indicated values of S. 



Figure 14: Models with different values of i5 are 
shown. As in the previous figures, for each value 
of L and Kq/L, only the best value is retained 
when the other parameters Vc and Va/y/Kg arc var- 
ied. The figure in the left panel displays the contour 
levels for < 40 for the indicated values of S. It 
is possible to scale the Kq/L values by a function 
f{6) to superimpose the contours corresponding to 
different values of S (see text). This is displayed in 
the right panel. 



Partial conclusion Several comments about 
these results can be made. First, the diffusion slope 
is constrained and the particular value S = 1/3 is 
excluded by our analysis (even if a less stringent cri- 
terion on x^ is adopted). We find that the power 
law index for the diffusion coefficient is restricted to 
the interval 0.45 - 0.85, the best x^ being 25.5 for 
5 = 0.70, leading to a '--^ 2.1. For any S in this in- 
terval, the good parameters in the Kq/L — L and 
Va/\^Ko — Vc planes can be straightforwardly de- 
duced from the corresponding values for S = 0.6 by a 
simple scaling law. We also exclude any model with- 
out a convective velocity or without reacceleration 
for any combination of the three other diffusion pa- 
rameters. These conclusions could get more stringent 
by new measurements in the whole energy spectrum 
for all nuclei, as all our results were obtained using 
the best data, which are rather scarce and more than 



5.3.2 Full exploration of the parameter space 

This study was then refined by adding four up- 
dates: i) a more conventional reacceleration term is 
used (i.e. similar to terms used by other authors); ii) 
effect of deviation from a pure power law of source 
spectra at low energy is inspected - see Sec. 2.4.3 , 
p. |ri| -; iii) different diffusion schemes (slab Alfvcn 
wave turbulence, isotropic fast magnetosonic turbu- 
lence and a mixture of the two case) are tested and, 
iv) above all, the source slope aj = a for all species 
is set as a free parameter. The net result, is that 
all the above trends are confirmed. As a bonus, one 
recovers that for small diffusion slope, convection is 
unfavored, as found in p7|. 



Summary of the trends In Fig. |6| we show the 
preferred values of the three diffusion parameters Kq, 
Vc and Va, for each best x^ in the S — j plane, when 
L has been fixed to 6 kpc (the behavior does not 
particularly depend on L). The two upper panels 
show that the evolution of a docs not affect Kq. On 
the other hand, we clearly see the (anti)corrclation 
between the two parameters Kq and S entering the 
diffusion coefficient formula, because they need to 
give about the same normalization at high energy 
{Kq X £^fhrcsh ~ ^^^) ■ Almost the same numbers are 
obtained for the pure power law and for t he mod ified 
source spectra (these were defined in Sec. 2.4.3 ). Kq 
spans the range 0.003 to 0.1 kpc^ Myr^^. The mid- 
dle panels show the values for the convective velocity. 
Only very few configurations include Vc = km s~^, 
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Figure 16: From top to bottom: for each best in the plane 5 
log(i4ro), Vc and Va are plotted for both source spectrum types. 



7 (L = 6 kpc), the corresponding values of 



tra. Increasing 7 and 6 at the same time makes Vc 
change its trend. Finally, the Alfvcn velocity Va dou- 
bles from S = 1.0 to 0.3, whereas it is almost un- 
changed by a variation in the parameter 7 (or equiv- 
alently a). B/C and sub-Fe ratios are displayed for 
several configurations in Fig.p7|. 

The three parameters Kq, Vc and Va behave very 
similarly with respect to a change in the source spec- 
trum from "pure power law" to "modified" . It can 
be explained as the influence on the primary and 
secondary fluxes which can be factored out if energy 
changes are discarded (their effect is actually small 
on the derived parameters). The currently available 
data on B/C do not allow to discriminate clearly 
between these two shapes for the acceleration spec- 
trum. This could be achieved by means of better 
data not only for B/C but also for primary nuclei 
(Donato et ai, in preparation). We finally notice 
that pursuing the analysis with the sub-Fe/Fe ra- 
tio does not allow to go farther for reasons invoked 
above. 



Conclusions Thus, forgetting for a while some of 
our theoretical a priori about the diffusion power 
spectrum, a new picture of cosmic ray propagation 
seems to emerge, in which high values for the diffu- 
sion coefficient spectral index (5 > 0.6 — 0.7) ^ and 
source spectral indices a ~ 2.0 are favored. The 
conclusions of the full analysis can be summarized 
as follows: (i) the values 6 ~ 0.7 — 0.9 and a ^ 2.0 
are preferred; (ii) this preference holds whatever the 
specific form of the spectrum at low energy; the nu- 
merical values of the other parameters are also only 
slightly modified by this low energy dependence even 
though deviation from a power-law at low energy is 
preferred. The study of fluxes should give a more 
definite answer; (iii) Kq scales logarithmically with 6 
and models with small halos tend to one-dimensional 
models with a simple relation between the surface 
mass density ^, Kq, L and Vc', (iv) several existing 
models are compared and the qualitative and quanti- 



*It is noticeable that Strong et al. that used and preferred 
S = 1/3 for the five nast vears took S = 0.42 — 0.52 in their 
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Figure 17: The B/C and sub-Fe/Fe spectra (modulated at cj) = 500 MV) for several sets of para meter s (giving 
the be st fit to B/C for these value s) are displayed, along wit h experimental data from heao-3 |14C| , balloon 
flights |141| , HEAT on ulysses |153| , hkh on isee-3 (Leskc |ll4| ) and voyager |154[| . Notice that ACE data [ll55[ 
correspond to a modulation parameter (f> « 750 MV. 



tative differences between them are studied and par- 
tially explained (the reader is referred to j24j for a 
deeper discussion). 



work has to deal with inspecting carefully the conse- 
quences of this degeneracy, and if possible, to break 
it. 



6 Diffusion parameters applied to 
"standard" CR physics 

Armed with our set of parameters derived from the 
B/C ratio, it seems reasonable to suppose that other 
charged nuclei follow the same propagation history; 
this is already the case for sub-Fe/Fe nuclei as shown 
in Fig. |l^. We therefore use the same diffusion 
model, with these sets of parameters, to study the 
propagation of all species. In this section wc focus on 
the subjects pertaining to the astrophysics field. Wc 
start with secondary p (Sec. 6.1) that is a further 
cross-check of the validity of the propagation param- 
eter derived. This nucleus is also particularly impor- 
tant to possibly constrain exotic primary sources in 
the diffusive halo |139, 152] - this will be presented 
in Sec. |^ -. We devote Sec. 3.2 to the question of 
the spatial origin of all these stable nuclei |^ . Then 
we focus on the study of radioactive species ||5^ in 
Sec. ^.3| , where some clues about the local properties 
of the interstellar medium will be given. Finally, in 



Sec. oA we study some consequences of our model 
on the composition at higher energies |6^ . 

It should be reminded that a strong degeneracy be- 
tween the propagation parameters was pointed out 



6.1 Secondary antiprotons and antideuterons 

The secondary antiprotons and antideuterons are 
yielded by the spallation of cosmic ray nuclei over 
the interstellar medium, mostly protons and helium 
whose spectra $p and <1>hc arc by now well deter- 
mined as discussed in Sec. 5.2.S . The production 
cross sections have been discussed in Sec. 2.4.6. 



6.1.1 Secondary antiprotons 

Unlike the secondary nuclei that are produced at 
fixed energy per nucleon, the source term for sec- 
ondary antiprotons is obtained from the convolution 
of the energy spectra of the incident cosmic rays with 
the relevant differential production cross sections. In 
the case of impinging protons, this leads to 



r + oo 

qf%T,Ep)=AT: dE^%{r,E,,) 

X <i "^H + n-Ho ,^ } (Ep Ep) 



dEp 



dEp 



A similar term arises from the cosmic ray helium 
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Tertiary production Once they have been cre- 
ated, antiprotons may interact with the intersteUar 
material in three different ways. First, they may un- 
dergo elastic scatterings on galactic hydrogen. The 
cross section for that reaction has been shown to 



peak in the forward direction [ |156[ so that the cor- 
responding antiproton energy loss is negligible. An- 
tiprotons are not perturbed by these elastic scatter- 
ings as they survive them while their energy does 
not change. Then, they may also annihilate on in- 
terstellar protons. This process dominates at low 
energy, and its cross section a^^^ is given in |157|. 
Last but not least, antiprotons may survive inelas- 
tic scatterings where the target proton is excited to 
a resonance. Antiprotons do not annihilate in that 
case but lose a significant amount of their kinetic 
energy. The cross section for these inelastic yet non- 
annihilating interactions is 



tPP 



non— ann 



(38) 



The energy distribution of antiprotons that have un- 
dergone such reactions has not been measured. It 
may be assumed to be similar to the proton energy 
distribution after p-p inelastic scattering. An im- 
pinging antiproton with kinetic energy Tp has then 
a differential probability 



dNp 



1 



(39) 



to end up with the final energy Ep, hence the differ- 
ential cross section 



dafi 



pH- 
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dEr, 
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(40) 



The corresponding source term for these so-called 
tertiary antiprotons may be expressed as 



qp''{r,E) = -47mHcr: 
-f 47rnH 



pp 

non— ann 



{E)%{r,E) (41) 
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The integral over the antiproton energy E of q^^' {E) 
vanishes. This mechanism does not actually create 
new antiprotons. It merely redistributes them to- 
wards lower energies and tends therefore to flatten 
their spectrum. Notice in that respect that the sec- 
ondary antiproton spectrum that results from the 
interaction of cosmic ray protons impinging on in- 
terstellar helium is already fairly flat below a few 
GeV. Since it contributes a small fraction to the fi- 
nal result, the effect under scruti ny h ere may not 
be as large as previously thought |158|. Helium is 
taken into account by replacing the hydrogen density 
in relation (^ij) by a geometrically-inspired factor of 
nH + 42/3nHc- 

Conclusions First, the values of all the inputs be- 
ing either extracted from the analysis of nuclei (dif- 



(proton and helium fluxes), all the cosmic antipro- 
ton fluxes naturally coming out of the calculation are 
completely contained within the experimental error 
bars of BESS data (see Fig. 18). 



> 
o 




Figure 18: Solid line shows the total top-of- 
atmosphere (TOA) secondary antiproton spectrum 
for the reference set of diffusion parameters (see text 
for details). Dashed lines are the contributions to 
this total flux from various nuclear reactions (from 
top to bottom: p-p, p-He, He-p and He-He). In data 
points, circles correspond to the combined BESS 1995 
and 1997 data ]l59C , squares to the 1998 ones [ pO | 
and stars to caprice data |l6l|| . 



The major uncertainties come from nuclear physics 
and are already comparable to observational error 
bars. In Fig. ^ the two dotted lines feature the 
uncertainties related to nuclear physics. The upper 
curve is obtained with the set of maximal p-He, He- 
p, He-He cross sections while increasing the p-p cross 
section - as given by the Tan & Ng parameterization 
- by a generous 10%. Similarly, the lower 



157| 



curve is obtained with the minimal values for those 
cross sections while decreasing the p-p cross section 
by 10%. Indeed, such a variation for p-p has been 
included for the sake of completeness even if it mod- 
ifies the antiproton spectrum only by a few pcrcents. 
The uncertainties which the upper and lower dot- 
ted curves exhibit are of the order of 22-25 % over 
the energy range 0.1-100 GeV. Another source of un- 
certainty comes from the fact that the propagation 
parameters are degenerate, even if they are severely 
constrained by the analysis of B/C experimental re- 
sults. This induces an error which is also displayed 
in Fig. ^ (solid lines). 

As antiproton spectrum measurements should bet- 
ter in the near future, antiproton studies could be 
limited by nuclear indeterminacies. Further work 
and especially new measurements of antiproton pro- 
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free from the effects of Solar modulation. This re- 
sult takes in to a ccount geomagnetic suppression as 
discussed in 162|. 



T- (GeV) 

Figure f 9: Solid lines show the envelope of the TO A 
antiproton spectra generated with the sets of dif- 
fusion parameters consistent with B/C. The dotted 
lines give an estimation of the uncertainty due to the 
indeterminacy of the nuclear parameters. 



est. 



6.1.2 Secondary antideuterons 

The production of secondary antideuterons also pro- 
ceeds through the collision of cosmic ray high-energy 
particles - mostly protons and helium nuclei - on the 
atoms of the interstellar medium. The source term 
q^'^ (r, Ejj ) for secondary antideuterons can then be 
readily obtained, using the coalescence scheme de- 
scribed in Sec. 2.4.7| . 

The result is that the spectrum of secondary an- 
tideuterons sharply drops below a few GeV/nuc (see 
Fig. ^ p. ^3|), mostly because in the galactic frame, 
the production threshold is 17 rrip. Energy losses 
tend to shift slightly the antidcuteron spectrum to- 
wards lower energies with the effect of replenish- 
ing the low-energy tail with the more abundant 
species which, initially, had higher energies. That 
effect is nevertheless mild. In particular, there is no 
tertiary contribution q^^ [r, Eq). An antideuteron 
that undergoes an inelastic scattering is likely to 
be destroyed unlike antiprotons for which inelastic 
yet non-annihilating interactions are possible (see 
Sec. |Tl]). 

As a consequence, the flux of secondary an- 
tideuterons is quite suppressed with respect to an- 
tiprotons. However, a preliminary study, which does 
not take advantage of the full set of propagation pa- 
rameters allowed by B/C shows that a dozen events 
may be collected by the AMS collaboration during its 
space station stage, in the energy range extending up 
to 100 GeV/nuc. These antideuterons correspond 



6.2 Spatial origin 

The question of the spatial origin of cosmic rays 
in our stationary diffusion model was addressed in 
Sec. ^. It was shown that the origin is fully de- 
scribed by three parameters L, Twind and rgpai (see 
also Sec. 4.3), which depend on the values of the 



propagation parameters. Having determined a defi- 
nite range of realistic parameters in Sec. |^, we can go 
one step further in characterizing the origin of cosmic 
rays. 



In Fig. 20 



Realistic values for rgpai and r„ 
are plotted Xw(<5, i) = 2i/r„i„d and Xspai(^,i) = 
i/rgpai for several energies and species. Actually, 
Xw{5,L) is a function of rigidity; it is more clever 
to use the latter choice instead of kinetic energy per 
nucleus, since Xw(^, i) then does not depend on the 
species. The left panel of Fig. |2^ displays Xw('^, L) for 
three rigidities: 1 GV, 10 GV and 100 GV. The quan- 
tity Xw {5, L) is an indicator of the competing role of 
convection and diffusion to keep the remote cosmic 
rays to reach Earth. The figure shows that above sev- 
eral tens of GV, diffusion has the main role, whereas 
convection may dominate below this value, at least 
for large values of 5. Different values of the halo 
size L yield the same conclusions. The right panel 
of Fig. |2^ displays Xspai('^, L) which is related to the 
spallation efficiency. It appears that at low energy, 
heavy nuclei such as Fe are preferentially destructed 
rather than swept out by the convective wind or by 
the escape through the boundaries. 



Origin of primaries and secondaries With the 
knowledge of rwind, ''spai and L at 1 GeV/nuc for all 
species one can now answer the question of the origin 
of CR: results for p, CNO and Fe nuclei are shown in 
Fig. ^ The peculiar behavior of secondaries is fully 
explained in [ p5| . 

To summarize the situation, most cosmic rays de- 
tected on Earth were emitted from a limited zone 
of the disc. This is even more true for heavy species 
such as Fe, which are very sensitive to spallations and 
thus are unlikely to travel long distances. This im- 
plies that the information on the diffusive processes 
inferred from the study of the ratio sub-Fe/Fc are 
only valid on a very local region, all the more that 
the diffusion slope 5 is large. Even if the propa- 
gation conditions were very different outside of this 
region, the observations made in the Solar neighbor- 
hood would almost not be affected, pointing prefer- 
entially once more towards large 5 values. Otherwise, 
as several anomalies in some CR radiations indicate, 
CR fluxes could be variable with location, making all 
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Figure 20: Left panel: Xw(<5, 7?.) as a function of the 
diffusion spectral index 6 for different rigidities TZ (it 
does not depend on species); from top to bottom, 
n = 100 GV, 7^ = 10 GV and 7^ = 1 GV. The 
parameter Xw, as well as Xspai, is not very sensitive to 
the halo size L. Right panel: XspaK*^? ^) as a function 
of S for 7^ = 100 GV (upper curves) and 7^ = 1 Gy 
(lower curves) for four species: p ((Tinoi ~ 40 mb), D 
(cTinei ~ 100 mb), B-CNO (crinci 250 mb) and Fe 
(CTinci ^ 700 mb). For the latter species we plotted 
the same three L values as in left panel. 



primaries (p, CNO, Fe) secondorles (B, sub — Fe) 
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Figure 21: Contours from where 99% of primaries 
(left panels) or secondaries (right panels) come, for 
two halo sizes L = 10 kpc (upper panels) or L = 
6 kpc (lower panels). From external to internal lines, 
and for each nucleus, these contours correspond to 
the diffusion slope values 5 — 0.35 — 0.60 — 0.85. 

6.3 The radioactive species 

The radioactive cosmic ray nuclei deserve a specific 
treatment. From their creation, they diffuse on a 



this expression, not only the diffusion coefhcient i^T, 
but also the lifetime 7T0, depend on energy, due to 
the relativistic time stretch. The following table 
gives some values of this distance for three species 
(typical values Kq = 0.033 kpc^ Myr^^ and S = 0.6 
have been assumed in Tab. H). These nuclei are 





TO (Myr) 


1 GeV/nuc 


10 GeV/nuc 


i^Be 


2.17 


220 pc 


950 pc 


26 Al 


1.31 


110 pc 


470 pc 


36 Q 


0.443 


56 pc 


250 pc 



Table 6: Rest frame lifetimes and corresponding val- 
ues of Irad for scvcral nuclei at selected energies. 

therefore very sensitive to the presence of the local 
bubble which has an extent of about 50 — 200 pc 
(see Sec. |Z|). As regards the diffusion process it- 
self, as described by the coefficient K{E), radio and 
7-ray observations, which can test in situ the spec- 
trum and density of cosmic rays 1 163 , indicate that 
it is not affected by the presence of the bubble, i.e. 
diffusion is homogeneous. 

The bubble has nevertheless an effect on the prop- 
agation. First it leads to a decrease in the spallation 
source term of the radioactive species. Second, it 
also leads to a local decrease of destructive spalla- 
tions. Third, as there is less interstellar matter to 
interact with, the energy losses are also lowered. Be- 
cause the typical propagation time of radioactive nu- 
clei is short, the energy redistributions are negligible, 
as well as destruction, and the first effect is domi- 
nant. As mentioned above, the bubble is modelled 
as a hole, in the thin disc approximation (see Fig. |ll], 
p. ^6|) . The radius of this hole is considered as an un- 
known parameter in the analysis. Fig. ^ shows, for 
''hole = 200 pc, how the fluxes both for stable and ra- 
dioactive species are affected in the neighborhood of 
the hole. The stables remain grossly unaffected by 
its presence, whereas radioactive arc strongly sup- 
pressed for the reason just mentioned. 

Effect and size of the hole The major result we 
found is that at the center of the bubble, the 
radioactive fluxes are decreased by a factor which 
can be approximated as 



jVholc 



;g OC exp(-rholc/'rad) • 



When all the sets of diffusion parameters allowed by 
the B/C data are used to compute the ^"Be/^Be, 
•^^Cl/Cl ratios, we find that each of the radioactive 
nuclei independently points towards a bubble of ra- 
dius < 100 pc, in relatively good agreement with 
direct observations. If these nuclei are considered 
simultaneously, only models with a bubble radius 
'"hole ~ 60 — 100 pc are consistent with the data. 
In particular, the standard case rhoic = pc is dis- 
favored. This is shown in Fig. |23l which is a pro- 
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Figure 23: Representation of the models compatible with B/C plus both "'^"Bc/^Be and '^^Cl/Cl ACE 3-cr (open 
circles) and l-cr (filled circles). Left panel displays homogeneous models (rhoic = 0) in the plane L — S. Right 
panel displays inhomogeneous models (fhoie > 0) in the plane rhoie — <5. 
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edge of :he hole (r„<,|. = 200 pc) 
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less clear , and it is suspected that the data (nu- 
clear or astrophysical) on which they rely should not 
be trusted. 



Whole energy spectrum As ^"Bc/^Be is prob- 
ably the best measured radioactive ratio (see 
Sec. 5.2.2), it is interesting to see how ratios de- 



pend on energy, given the constraint that parame- 
ters must fit both B/C and ^"Be/^Bc data. This is 
shown in Fig. |2^ for the homogeneous model and a 
hole Thole = 80 pc. These curves show that it would 
be very important to have accurate measurements 
of the energy spectra of radioactive species such as 
lOBe, 36C1 and ^^Al. 



r(kpc) 



Figure 22: Radial distribution of a radioactive 
species in the disc, across the hole, for rhoio = 200 pc. 
Numerical values were taken for ^^Al (radioactive) 
and ^^Si (stable). The distribution of radioactive is 
very sensitive to the presence of the hole. However 
this effect is local and vanishes at several rhoic- On 
the other side, the distribution of stable species is 
not much affected. 



lOBe/^Be and ^^Cl/Cl on the L - J plane (left panel, 
no hole) or rhoio — ^ plane (right panel, hole J'hoio)- 
When the ratio ^''A1/^''A1 is considered along with 



Conclusions Most studies use the radioactive nu- 
clei to constrain the halo size L (sec discussion in 
Sec. 2.4.2). The meaning of the results are very much 
dependent on the the treatment of the LISM. In 
particular, the presence of the local bubble leads to 
an exponential attenuation of the radioactive fluxes. 
This is of special importance for short-lived nuclei 
such as ^^C, which is attenuated by a factor ^ 1, 
unless a local source is present. This may explain 
the fact that ACE did not detect any such nucleus. 
Finally, if one takes rhoic as an input fixed by obser- 
vations, the radioactive progenitor /parent may give 
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Figure 24: Envelopes of the spectra obtained with all the models compatible with ^°Be/^Be ACE 1-a for 
the three radioactive species. Solid lines are for homogeneous models (rhoie = pc) and dashed lines are for 
inhomogeneous models (rhoio — 80 pc). Data are from ACE (circles), ulysses (crosses), voyager (filled squares) 
and ISEE (see Sec. 5.2.2). 



6.4 PeV fluxes 

Recent measurements of the cosmic ray average 
logarithmic mass and all-particle spectrum around 
10^^ eV |164, 165| gave new clues to understand the 



origin of the cosmic rays and in particular the puzzle 
of the knee in the energy spectrum. The highest en- 
ergy particles are almost certainly extragalactic. A 
similar origin is not excluded near the knee, but it is 
difficult to account for the observed softening of the 
spectrum in this region. As a consequence, the inter- 
mediate region between 10^^ and 10^^ eV should be 
analyzed in terms of the same physical mechanisms 
than lower energy particles. 

At these energies, fluxes are too low to be di- 
rectly measured (~ m~^ sr~^ yr~^). With present 
techniques, only two quantities can be extracted in 
large ground array detectors (e.g. |166 ); namely all- 
particle spectrum and (In A) . They are given by lin- 
ear combination of the individual fluxes with differ- 
ent weights; 



all _ 



, and (lnA)^ % ' .(42) 



As the experimental data are given in total energy 
rather than in kinetic energy per nucleus, we adopt 
the same presentation for the results above 100 GeV. 



6.4.1 Separation of key ingredients 



some basic ingredients that drive the mass evolution, 
but not all of them. In particular, the geometrical 
effects cannot be considered. Trends are also more 
easily understood with only one light and one heavy 
nucleus (e.g. p and Fc). If energy gains and losses 
arc discarded, the Leaky Box equation for primary 
species, such as p and Fe, reads 



+ qlQ^{E)-rN^ =0 



(43) 



The proportion of light and heavy species evolves 
with energy because of two effects. The source spec- 
tra ql^Q^{E) may differ, and the spallative destruc- 
tion rate F^ are definitely not the same. The escape 
time Tosc = tqRt^ is supposed to be independent of 
species. 



Source spectrum effect The first effect is stud- 
ied by neglecting the spallation term in Eq. (p3|). 
We suppose that the two species p and Fe are in- 
jected with different slopes related by apc — Q^p = 
0.1 [ [167| . If we start from the experimental value 
(Fc/p)ioo Gcv = 1/20 (e.g. |l6|), Eq. (g) gives the 



following evolution, independently of 5, 

(In A){E = 100 GeV/ 10 TeV/ 1 PeV/ oo) 
~ 0.19/ 0.30/ 0.50/ 4.03. 



We first estimate the evolution of (In A) in the Leaky 



Escape plus selective destruction If the spalla- 
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as 



{lnA)iE)^HAF,) X 



l + (Mc)x 



p 
-, -1 



{aoiE/Z^y 



(44) 



where the cross sections are dp""' 



30 mb and 

710 mb). The quantity Tesc has been con- 
verted in mb, with a typical value ctq = 47 mb 
and S = 0.61 We are not interested here in 

the effect of injection spectra and the slopes have 
now been assumed to be the same for all species 
(ap = apo = a=2.2). With (In ^) (100 GeV) = 0.19, 
this gives 

{lnA){E = 100 GeV/ 10 TeV/ 1 PcV / oo) 
- 0.19/ 0.83/ 1.05/ 1.08. 



Geometrical effects are subdominant Fig. 25 

shows spallations effects and geometrical effects. The 
first strong conclusion is that spallation effects af- 
fect dramatically the composition of cosmic rays. 
Furthermore, at sufficiently high energy, as expected 
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Geometrical source effect As emphasized 
above, the study of the influence of geometry is not 
straightforward in Leaky Box models. From the 
previous discussion in the framework of diffusive 
propagation (see Sec. ||), we know that cosmic 
rays detected on Earth come from an average 
distance (r){E) which is larger at higher energy. 
As a metallicity gradient is present, at the level of 

0.05 dex kpc-i, i.e. [Fe/H] = -0.05 (r - Rq) = 

log (gpc/gp)''*^ - log(Fe/H)0 (see Sec. |2.4.3D , the 
nuclei with different energies were emitted from 
sources having different properties. This leads to 
(the two previous effects are ignored) 



{\nA){E) = 



ln(A 



Fc) 



1 + (9p/9Fo) 



eft 



where 



cff 



An upper bound is obtained for this effect by as- 
suming that (r)(100 GeV) = Rq and (r)(oo) = 
(galactic center), which yields 

(In A) {E = 100 GeV/ oo) - 0.19 / 0.54 . (45) 

Even in this optimal case, this effect is seen to be 
negligible compared to the two others (see below). 



6.4.2 Results from our propagation model 

The quantity (In A) can also be computed with our 
propagation model, considering the different sets of 
propagation parameters derived from the B/C anal- 



Figure 25: Average logarithmic ma ss (In A) for the 
radial distributions depicted in Sec. ^.4.3| , p. [lO| ((5 = 
0.6 and L = 10 kpc for all models). 

with the Leaky Box toy model, the asymptotic 
regime is reached and propagation ceases to affect 
(in A). The geometrical effects (source distribution 
plus metallicity gradient) are less important than the 
others. They induce a change of at most 5% in the 
results. This would be even more true for lower L. 
For L = 3 kpc, these geometrical effects are com- 
pletely negfigible. 



Evolution above PeV energies The evolution 
of (In A) obtained with different sets of parameters 
compatible with B/C is presented in Fig. |2^. The 
right panel focuses on the source effect and shows 
that Op > a is preferred. As suggested in our toy 
model, only source effects enable evolution of (In A) 
above 1 PeV. It can also be seen that large diffu- 
sion slopes S are preferred, as was already hinted in 
Sec. 



5.3 



Inclusion of the knee An important experimen- 
tal fact, referred to as the knee, has not be taken 
into account in the analysis presented above. It is 
a drop of the slope of A7 — 0.4 in the observed 
spectrum of cosmic rays at an energy of a few PeV. 
This may be due to a change in the source spec- 
tra or in the diffusion coefficient. We now want 
to study the effect of the knee on the evolution of 
(In A). To this aim, we propose to model the knee 
in a very naive way, by a break in the source spec- 
tral indices at an energy -Eknoo that depends on the 
species (notice that a change in d iffusive regime, e.g. 
as recently re- inspected in |16£], is also possible). 
We consider two break at a given rigidity 

-£'knee = Z X i PeV and a break at given energy per 
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Influence of 5 and a (L=10 kpc) 




Figure 26: For different sets of parameters compat- 
ible with the B/C ratio, the average mass {lnA){E) 
normahzed to 100 GeV observations is displayed as a 
function of energy. Upper panel: the diffusion slope 
d is varied keeping aj = a for all species (thin lines, 
no geometrical effects; thick lines, geometrical effects 
are included). Lower panel: a different source slope 
is considered for p {S is set to 0.6). Stars show ex- 
perimental data [164| (see [p8| for details). 



mass {In A) is shown in Fig. (the all particle flux 
can be found in [^). Several features are vis- 
ible. They are first generated by protons changing 
their spectral index, followed by others: the helium 
component, the CNO group, and Fe. At high energy, 
all species have the same post-knee spectral index. 
The transition is even more gradual when secondary 
species (upper curve) are taken into account (they 
also change the normalization). The choice of the 
break (in rigidity or in energy per nucleon) affects 
mostly the all-particle spectrum. The conclusion at 
this stage is that the scatter in the data is too large 
to provide a clear test of our theoretical predictions 
(see Fig. 



7 Propagation models applied to as- 
troparticle physics 

Connections between high-energy cosmic ray physics 
and cosmology are explored in these three sub- 




E [GeV] 

Figure 27: Average logarithmic mass for models in- 
cluding a break either in rigidity (i? = 4 PeV) or in 
total energy per nucleus (E/nuc = 4 PeV). Solid lines 
correspond to primaries only, whereas dotted lines 
correspond to primaries plus secondaries. For illus- 
trative purpose, some data from ground arrays have 
been displayed: triangl es ar e Casablanca's data 
with Monte-Carlo hdpm |17C|, empty circles are from 
KASCADE collaboration with QGS jet simulation [ll66[ 
(empty squares are from jacee direct experiment); 
RUNJOB data plus some jacee data are lower than 
1.8 and do not appear on the graph, see 164|). 



lustrate the above-mentioned ties by a discu ssion (i) 
on the limits on the cold invisible gas [171|, (ii) on 



the baryon asymmetry in the universe (see e.g. [172|) 
and finally (iii) on the presence of putative neutrali- 
nos that could make up the astronomical dark mat- 
ter 1 173 1 or primordial black holes that could have 
been formed in the early Universe. Then, in Sec. 7.2 
we explore the question of the spatial origin of these 
exotic nuclei. Finally, Sec. 7.3 and Sec. |7.4| respec- 



!H1|, 



139| 



tively deal with SUSY [ll52|, |l6g and PBH 
signatures in antiprotons and antideuterons. 



7.1 Introduction 

Ties are strong between high-energy cosmic ray 
physics and cosmology. In this section, we explore 
a few examples that illustrate these connections and 
show that a deep understanding on how cosmic rays 
propagate turns out to be crucial. 

A long standing problem in astronomy lies in the 
existence of large amounts of unseen material whose 
gravitational effects are nevertheless large. To com- 
mence, the Milky Way is surrounded by an extended 
halo that induces a flat rotation curve in its plane 
\174\. This trend has also been observed in many 
spiral systems. On larger scales, the presence of 
dark matter inside galac tic c lusters has been no- 
ticed since many decades 175 . Finally, the observa- 



tions of the Cosmic Microwave Background (CMB) 
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nation of the relation between the distance of himi- 



nosity and the redshift of type la supernovae ||l77| , 
or with the large scale structure (LS S) information 
from galaxy and cluster surveys [ [178| , give indepen- 
dent evidence for a cosmological dark matter density 
in the range ^Icdmh^ ~ 0.13 ±0.05 |176|; to be com- 
pared to a baryon density of fib — 0.019 ± 0.002 



as indicated by nucleosynthesis [179 and the relative 
heights of the first acoustic peaks in the CMB data, 
= 0.022 ±0.003 [|l8§. The nature of the astro- 
nomical dark matter is still unresolved and most of it 
could be made of non-baryonic species or primordial 
black holes. But observations also point towards the 
presence of dark baryons insofar as the amount of 
luminous material is significantly smaller than fib- 



Limits on the Dark Gas inside the Milky Way 



As a dark matter of fact. Dc Paolis et al. |181 have 
outlined a scenario in which dark clusters of com- 
pact objects pervade the halo of the Milky Way to- 
gether with clouds of molecular hydrogen, at dis- 
tances larger than 10 to 20 kpc. Pfenniger et al. [ |l82| 
have also suggested that a flat rotation curve beyond 
the Solar circle could be explained by a thin disc of 
cold molecular hydrogen, widening at large distances 
from the center. If unseen gas was present in the 
galactic halo or in the ridge of the Milky Way, it 
would be impacted by cosmic rays originating from 
the disc. This would lead to a strong 7-ray signal 
showing up as a new component in the galactic dif- 
fuse radiation. By using the C OSB results at high 
galactic latitude, Gilmore |183| has already inferred 
a conservative limit ~ 15 % on the amount of hidden 
gas in the halo. That analysis can be significantly re- 
fined because a key ingredient - the propagation of 
cosmic rays - is now better understood. A reliable 
calculation of their density as a function of galac- 
tocentric radius r and height z above the plane is 
definitely possible, as discussed in the previous sec- 
tions. It allows for the determination of the gamma 
ray hydrogen emissivity Ir per hydrogen atom as a 
function of location 



lH{E^,r,z) 



W(r0,O) 



Ih{E^,Q) 



where the spectral index of cosmic ray protons - the 
dominant contribution - is found to remain roughly 
constant throughout the Galaxy. In the local range, 
the emissivity is 



lH{E~f > 100 MeV,0) 



(1.84 ± 0.10) X 10"^*' photons H"^ s"^ sr^^ 

The spallation of cosmic rays with the gas potentially 
concealed in the Milky Way produces an extra 7- 
ray diffuse emission whose flux obtained from the 
convolution along the line of sight of the density n^'^^ 
of dark gas with the local 7-ray emissivity Ih is 
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Figure 28: A Mestel's disc generates the flat galactic 
rotation curve. The central part of that disc, up to 
the radius rmin, does not contain any cold gas so far 
undetected in CO. In the left panel, the inner gas 
boundary rmin is varied for various values of the disc 
thickness Odisc- Conversely, in the right panel, the 
bound ?/max is presented as a function of Odisc [171|. 



Observations collected by the egret instrument on 
board CGRO have been compared to the 7-ray diffuse 
radiation that may be modeled with both HI and II2 
contributions. An isotropic component is found in 
the residuals of the signal. Its integral flux between 
500 MeV and 1 GeV is found to be 2.5 x 10"^ photons 
cm^^ s^^ sr^^ in the direction of the Perseus arm 
(b = 0° , 1 = 115°). By requiring that does 
not exceed that value, a constraint may be derived on 
the amount of cold molecular clouds - so far untraced 
by CO - which the Milky Way ridge could conceal. A 
Mestel's disc with local surface mass density E(0) ~ 
200 Mq pc~^ is assumed. The spatial density of that 
distribution is given by 

-E(0) 



(47) 



where Odisc denotes the disc thickness. Fig. ^ fea- 
tures the upper limit Tymax on the fraction of cold 
gas which the outer parts of the disc (r > rmin) may 
contain. In particular, for Odisc = 400 pc, the limit 
relaxes to 1 only if rmin exceeds 19.5 kpc. We there- 
fore conclude that molecular hydrogen, sufRciently 
cold to have escaped detection and so far untraced 
in CO, should actually be sparse up to ~ 20 kpc from 
the galactic center and that it cannot account for the 
rotation curve of the Milky Way. 



Neutralino as a Dark Matter Candidate, PBH 



7.2 Spatial origin of the PBH or SUSY exotic primaries 
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candidates to the astronomical missing mass is a neu- 
tral weakly interacting particle. Such a species is pre- 
dicted in particular by supersymmetry, a theory that 
is actively tested at accelerators. There are indirect 
clues for supersymmetry, noticeably the existence of 
a single high-energy unification scale for supersym- 
metric grand-unified models alone. It is conceivable 
therefore that most of the dark matter in the halo 
of the Milky Way is made of such neutral particles. 
The relic abundance of these so-called neutralinos is 
relevant to cosmology. If present in the galactic halo, 
they should still annihilate mutually to yield, among 
a few other indirect signatures, a flux of antiprotons 



whose spectrum is discussed in Sec. 7.3 



The same line of argument applies if black holes 
make up a fraction of the galactic dark matter. Pri- 
mordial black holes could have formed in the early 
universe and could have concealed baryons that are 
not accounted for by primordial nucleosynthesis since 
they would not have participated in the nuclear fu- 
sion of helium. The actual value for fib could there- 
fore be much larger than ~ 5% and even reach llcdm) 
sparing us the need for non-baryonic particles. If the 
black holes are light, they could evaporate today and 
yield also antiprotons as is discussed in Sec. 7.4. 



The primary antiprotons produced by neutralino 
annihilations or black hole evaporations will get 
mixed with a conventional population of secondary 
antiprotons that are produced by the spallation of 
cosmic ray nuclei on the interstellar gas of the Milky 
Way ridge. It is crucial to investigate in some detail 
this later mechanism in order to ascertain which of 
the primary or secondary species is dominant. This 



has been done in Sec. 3.1 



Antinuclei in the Galaxy Another crucial prob- 
lem of cosmology is the existence of a possible asym- 
metry between matter and antimatter. The current 
wisdom is that we live in a matter dominated uni- 
verse and that antimatter islands seem to be ex- 
cluded. The ongoing annihilations that would take 
place at the frontiers would create an intense gamma- 
ray emission that is not seen. The mechanism which 
would be required to account for a separation of mat- 
ter and antimatter domains is not known. However, 
one should not subscribe to any dogmatic point of 
view. The amount of antimatter in cosmic rays is 
about to be measured with unequalled accuracy by 
the space station borne spectrometer of the AMS col- 
laboration [184|. One of the most exciting goals of 



the experiment is the possible detection of antinu- 
clei in the cosmic radiation. It is generally believed 
that the observation of a single anti-helium or anti- 
carbon would undoubtedly signal the presence of an- 
tistars because the conventional production of antin- 
uclei through spallatio n is negligible. This point has 
been examined in Sec. xl where we mostly concen- 
trated on secondary antideutcrons. Suffice it to say 
that the creation of any additional antinucleon dur- 



of - 10-5 - 10-4 ||. As featured_ by Fig. |§, the 
D/p ratio is - 3 x 10-^ while Ke/D - 7 x 10"^ for 
a momentum per nucleon in excess of 10 GeV. Sec- 
ondary antinuclei with atomic number A > 3 are so 
much suppressed that they are out of reach of the 
near future instruments. Another possibility would 
be that antiglobular clusters lies in the halo of the 
Galaxy and could make a substantial contribution 
to the anti- helium flux |185]. There is more hope 



for antideutcrons. Their primary production will be 



discussed in Sec. 7.3 for supersymmetric neutralinos 



and briefly mentioned in Sec. 7.4 in the case of evap- 
orating black holes. 
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Figure 29: The fluxes of cosmic ray antiprotons and 
of anti-deuterium and anti-tritium nuclei, relative to 
the proton flux, are presented as a function of the 
momentum per nucleon | |9^ . To flt on the same di- 
agram, the curves have been scaled by a factor of 
10^ for anti-deuterium and of 10^ for anti-helium He. 
The doubling of curves corresponds to different fac- 
torization schemes. 



7.2 Spatial origin of the PBH or SUSY exotic 
primaries 



As for species created in the disc (see Sec. 6.2), one 
can infer the closed surfaces from which most of the 
cosmic rays detected on Earth ]l86t were emitted 
from. We are now confronted with a volume distri- 
bution of sources (see Sec. 2.4.4). As for standard 



sources in the disc, geometrical effects arc very im- 
portant (see Sec. |4.3| ). Apart from these geometrical 
restrictions, spallations and convectivc wind greatly 
affect the shape of these isodcnsity surfaces and, e.g., 
location of the particularly interesting contour de- 
fined by V(y{x)\0) = 99% (contour from where 99% 
of exotic cosmic ray originate). As a consequence, 
not only the characteristic size L is important to 
estimate the x = 99% surface, but also spallations 
and convection typical extension, i.e. parameters 
Tspai and Twind such as defined in Eqs. ( |3^ ) and (^5|). 
These parameters strongly depend on 6 (the diffu- 
sion coefficient slope) when realistic configurations 
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Realistic values of rgpai and r-wind for an- 
tideuterons Aiitideuteron signal seems to be the 
most promising species to consider. We choose the 
interstellar energy 1 GeV/nuc; nuclei that reach the 
detector are Solar modulated so that the final en- 
ergy corresponds roughly to 200 — 800 McV/nuc de- 
pending of the modulation parameter. This is the 
window where the exotic signal becomes interesting. 
Tab. 1^ summarizes the values of r^ind and rgpai at 
this energy for antideuterons for three halo sizes and 
three values of S. We notice that the situation is 





(kpc) 


S = 0.35 5 = 0.6 6 = 0.85 


L = 10 kpc 


^wind — 
^spal — 


OO 8.95 3.29 
24.19 8.69 4.02 


L = 2 kpc 


^wind — 
^spal — 


oo 2.42 0.95 
6.97 2.49 1.18 



Table 7: Twind and rgpai for two halo sizes L and three 
diffusion slopes S: these numbers are for 1 GeV/nuc 
(interstellar energy) antideuterons. 



very different for small or large S. For Kolmogorov 
power spectrum, only spallations act on the prop- 
agation and only weakly for this light nucleus: for 
large i5 - the value S = 0.85 is the one preferred in 

models are convec- 



our B/C analysis (see Sec. 3.3) 
tion/spallation dominated with spallation and con- 
vection acting at about the same footing. 
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Figure 30: Contours defining surfaces from where 
99% of exotic primaries come (no side boundaries). 
Upper panels: cut in the 2^ = kpc plane; lower 
panels: cut in the = kpc plane. Left panels 
correspond to L = 2 kpc and right panels to L = 
10 kpc. In each panel, we plot either the PBH case 
(solid lines) or the SUSY case (dotted lines). From 
external lines to internal lines correspond the values 
of the diffusion coefficient slope 8 = 0.35, i5 = 0.60, 
5 = 0.85. 



Results The production of particles by PBH is 
related to the dark matter profile whose weight is 
wp-Bnir, z) = f(r,z), whereas for supersymmetric 
annihilating particles, the production is rel ated t o 
the weight wsusy('", z) = f{r,z)^ (see Sec. 2.4.4 ). 
Weighting the elementary probability by this produc- 
tion rate leads to the contours displayed in Fig. |3^ 
for 7^cyi(V(a::)|0) = 99% (parameters arc taken from 
Tab. 0). For S ~ 0.35 (external contours), we re- 
cover basically contours that one would obtain in the 
high energy limit, i.e. contours driven by geometri- 
cal effects, pure diffusion (see also [Q). However, 
for larger S (internal contours), these contours shrink 
naturally and all surfaces are distorted towards the 
galactic center, where the maximum of production 
occurs. One has also to keep in mind that whatever 
these surfaces, most cosmic rays emitted from inside 
these volumes actually do not reach us. Escape pref- 
erentially occurs and it is obvious that the closer the 
source, a smaller part of the emitted flux escapes. 
The fraction that reaches us can be estimated to be 
about 0.01%-0.1% of the total emission in volume de- 
fined by a; = 99%, whatever S (it slightly depends on 
the halo size L for PBHs as well as for SUSY parti- 
cles). Another meaningful number is the fraction of 
particles created in the volume x ~ 99% with respect 
to the total number of particles created in the entire 
dark halo. This quantity is strongly dependent on 
the diffusive halo size and, for L = 10 kpc compared 



80% vs 20% (SUSY). This is a mere consequence of 
the different production terms. 



7.3 Primary antiprotons and antideuterons 
from supersymmetric sources 

The neutralinos that could be concealed in the halo 
of the Milky Way - and be responsible for the flat- 
ness of its rotation curve - should be steadily anni- 
hilating and produce antiprotons together with an- 
tideuterons. A key difference with respect to the 
production mechanisms that have already been dis- 
cussed lies in the fact that this new yet putative 
source of antinuclei is not confined to the galactic 
disc. It spreads all over the halo far above and be- 
neath the ridge. Its contribution to the local flux 
may nevertheless be easily derived with the formal- 



ism presented in Sec. 3.2.3. The propagation of pri 



mary species from the remote regions of the Milky 
Way neighborhood to the Earth has actually been 



treated in |152] 



SUSY model The neutralino naturally appears in 
the framework of the Minimal Supersy mmet ric ex- 
tension of the Standard Model (MSSM) ITs^ as the 
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the zino Z and the two higgsino H° and states 

X = ai7 + + agi/r + 04^2° ■ (48) 

This particle is neutral and interacts weakly. It gen- 
erally turns out to be the lightest supersymmetric 
state and is therefore stable if R-parity is conserved. 



It may be expressed as a sum - extending over the 
various quarks and gluons h as well as over the dif- 
ferent annihilation channels F - of the square of 
the antiproton differential multiplicity. That sum 
is weighted by the relevant branching ratios. The 
antineutron and antiproton differential distributions 
have been assumed to be identical. This readily leads 
to the source term for supersymmetric antideuterons 



Antiproton multiplicity The differential multi- 
plicity for antiproton production in a neutralino pair 
annihilation may be expressed as 



dNp 



E 

F,h 



B 



(F) 



dNl 
dEr, 



The annihilation proceeds - through the various 
final states F - towards the quark or the gluon h 
with the branching ratio . Quarks or gluons 
may be directly produced. They may alternatively 
result from the intermediate production of a Higgs or 
a gauge boson as well as of a top quark. Each quark 
or gluon h generates in turn a jet whose subsequent 
fragmentation and hadronization yields the antipro- 
ton energy spectrum dNp /dEp. The source term for 
supersymmetric antiprotons 



ql'^''^{r,Z,Ep) = (cTannW) 



dNp f p^{r,z) 
dEr, 



supplements the spallation contributions g^^"^ and 

Because the annihilation of a neutralino pair x~X 
occurs at rest with respect to the Galaxy, most of the 
antiprotons - and antideuterons for that matter - are 
produced at low energies. The resulting spectrum 
is fairly flat below a few GeV. This has important 
observational consequences insofar as the secondary 
antiproton spectrum is already fairly flat as shown in 
Fig. [l^, p. HJ. It may therefore prove difficult to dis- 
entangle a primary component from the secondary 
antiproton radiation. The primary flux nevertheless 
modifies the magnitude - if not the spectrum - of 
the cosmic ray antiproton radiation and could con- 
tribute significantly as shown in Tab. ^. This effect 
can be used to constrain the supersymmetric param- 
eter space. 



Supersymmetric D signal As regards the an- 
tideuteron production, the factorization-coalescence 
scheme discussed above leads to the antideuteron dif- 
ferential multiplicity 
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Figure 31: The IS flux of secondary antideuterons 
(heavier solid curve) decreases at low energy whereas 
the energy spectrum of the antideuterons from super- 
symmetric origin tends to flatten (see Sec. |7.3| ). The 
four cases of Tab. ^ are respectively featured by the 
solid (a), dotted (b), dashed (c) and dot-dashed (d) 
curves. Solar modulation has been taken at maxi- 
mum when the AMS observatory operates on board 
ISS. 

The four supersymmetric configurations of Tab. js] 
are presented in Fig. |l] together with the calculated 
secondary spectrum (see Sect.6T). The antideuteron 
spectra are fairly flat at low energy. They yield a few 
events for the AMS observatory on board ISS below 
an IS energy of 3 GeV/nuc - a region from which 
secondary antideuterons are absent -. 



Consequence for SUSY parameter space For 

each configuration of the whole supersymmetric pa- 
rameter space, the D flux has been integrated over 
that low-energy range. The resulting yield Nq which 
AMS may collect is presented as a function of the neu- 
tralino mass in the scatter plot of Fig. During 
the AMS shuttle mission, the Solar cycle was close to 
maximum. Most of the configurations are gaugino 
like (crosses) or mixed combinations of gaugino and 
higgsino states (dots). A significant portion of the 
parameter space is associated to a signal exceeding 
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case 








(0.24 GeV) 


(0.24 GcV/mic) 


$'5^^^ (0.24 GeV/nuc) 


^vrmax 
D 


a 


36.5 


96.9 


0.20 


1.2 X 10"^ 


1.0 X 10-'^ 


2.9 X 10-« 


0.6 


b 


61.2 


95.3 


0.13 


3.9 X 10-3 


3.5 X lO-'^ 


1.1 X lO-'^ 


2.9 


c 


90.4 


53.7 


0.03 


1.1 X 10-3 


1.8 X lO-'^ 


6.1 X 10-^ 


2.0 


d 


120 


98.9 


0.53 


2.9 X 10--* 


2.5 X 10-^ 


8.6 X 10-^ 


0.3 



Table 8: These four cases illustrate the richness of the supersymmetric parameter space. There is no obvious 
correlation between the antiproton and antideuteron Earth fluxes with the neutralino mass m-)^^. Case (c) 
is a gaugino-higgsino mixture and still yields signals comparable to those of case (a), yet a pure gaugino. 
Antideuteron fluxes are estimated at both Solar minimum and maximum, for a modulated energy of 0.24 
GeV/nuc. The last column features the corresponding number of D's which AMS on board ISS can collect below 
an IS energy of 3 GeV/nuc. 
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Figure 32: The supersymmetric D flux has been 
integrated over the range of IS energies extending 
from 0.1 up to 3 GeV/nuc. The resulting yield Nj^ 
of antideuterons which AMS on board ISS can collect 
is plotted as a function of the neutralino mass m^. 
Modulation has been considered at Solar maximum. 



few cases, AMS may even collect more than a dozen 
of low-energy D nuclei. Notice finally that numer- 
ical simulations have shown that ncutralinos - and 
more generally cold dark matter - should cluster in 



very dense and numerous clumps [188|. Because neu- 
tralino annihilations proceed through a two-body re- 
action, their antiproton and antideuteron signatures 
would be enhanced by a large factor that could even 
reach up a few hundred in the case of density pro- 
files d la Moore. If so, the entire supersymmetric 
constellation in Fig. ^ would be shifted upwards by 
at least two orders of magnitude and detection would 
become crystal clear for the most optimistic config- 
urations. Another important question is related to 
the dumpiness of the dark halo. If the dark mat- 
ter is concentrated into clumps, then the cosmic ray 
signals on Earth may be sensitive to the spatial dis- 
tribution of the clumps. This aspect can be studied 
with the formalism used to determine the spatial ori- 



in preparation). 



7.4 Primary antiprotons and antideuterons 
from Primordial Black Holes 

Very small black holes should have formed in the 
early Universe from initial density inhomogeneities 
(Hawking |189 ). They should now evaporate in- 
tensely through the Hawking mechanism (Hawking 
]19C| |) if their initial masses were around w 
5 X 10"'^'* g. Detecting such objects is a great challenge 
of modern physics and cosmology as it would both 
allow to give experimental grounds to the Hawking 
radiation which is one of the only tentative achieve- 
ment of semi-classical quantum gravity, and to probe 
the very small scales of the early Universe that re- 
mains totally inaccessible to other observations. Al- 
though the standard cosmological model of structure 
formation, assuming a pure scale-invariant Harisson- 
Zcldovitch power spectrum normalized to CMB am- 
plitudes, would lead to a very small amount of PBH 
dark-matter in the present Universe, several realistic 
infiationary scenarios and phase-transition phenom- 
ena can produce a significant amount of PBHs. 



7.4.1 Hawking evaporation and antiproton source 
term 

The Hawking black hole evaporation process can be 
intuitively understood as a quantum creation of par- 
ticles from the vacuum by an external field. The ba- 
sic characteristics can be easily seen through a sim- 
plified model, and the interested reader is referred 
to [191 for more details. 



Elementary energy spectrum The accurate 
emission process was derived by Hawking, using the 
usual quantum mechanical wave equation for a col- 
lapsing object with a post-collapse classical curved 
metric. He found that the emission spectrum for par- 
ticles of energy Q per unit of time t actually mimics 
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dom: 



dQdt 



cxp 



Q 



(49) 



where the effect of angular velocity and electric po- 
tential have been neglected since the black hole dis- 
charges and finishes its rotation much faster than it 
evaporates [192 . In this expression, k denotes the 
surface gravity, s is the spin of the emitted species 
and Tg is the absorption probability, given by 



47rcrs(Q,A/,/i) 2 



(50) 



where as is the absorption cross section computed 
numerically |193| and /x the rest mass of the emit- 
ted particle. It is also convenient to introduce the 
Hawking temperature, defined by 



T : 



he" 



16nkGM 



1013 g 
M 



GeV . 



(51) 



The expression ( p9[ ) may then be written simply as 
a function of Q/kT. 

Antiproton multiplicity As it was shown 
in 1 194 1, when the black hole temperature is greater 



than the quantum chromodynamics confinement 
scale Aqcd, quarks and gluons jets are emitted in- 
stead of composite hadrons. To evaluate the number 
of emitted antiprotons p , one therefore needs to per- 
form the following convolution: 



dM/dt = —a{M)/A'P (by a simple integration of 
the Hawking spectrum multiplied by the energy of 
the emitted quantum) where a{M) accounts for the 
available degrees of freedom at a given mass. In the 
assumption a{M) « const it leads to: 



d^n{M) _ 
dMdV^ ~ {iat + M3)2/3 ' dMdV^ 



Mi=M? 



(53) 



where the derivative is evaluated at = {3at + 
M3)1/3. The initial values d'^n,/dM,dVi were 
shown |196| to scale as for PBHs formed, as 

expected, in a radiation dominated Universe from a 
scale invariant power spectrum. The nowadays spec- 
trum is therefore mostly identical to the initial one 
above A/» = Sato ~ 5 x lO^^g and proportional to 
below. 

The last problem is that the spatial distribution 
of these PBHs is basically unknown. However, these 
objects should have formed in the very early stages 
of the history of the Universe: their initial mass is 
very close to the horizon mass at the formation epoch 
{Mh w Mpi X tfoim/^pi which means that the forma- 
tion time tform is of the order of 10^'^'^ s for M ~ A'h). 
Once formed, they only interact through gravitation, 
so that they behave as cold dark matter, and their 
spatial distribution should be similar. As a conse- 
quence, we use the same profile for the PBHs distri- 
bution as for the neutralinos, as given in Eq. ([l5|). 

7.4.2 Resulting upper limit and cosmological conse- 
quences 



dEdt 



E 



rs,(g,r) 



Q=E 



dgjpiQ^ 



where Uj is the number of degrees of freedom, E is 
the antiproton energy and dgjp{Q, E)/dE is the nor- 
malized differential fragmentation function, i.e. the 
number of antiprotons created with an energy be- 
tween E and E + dE by a parton jet of type j and en- 
ergy Q. The fragmentation functions have been eval- 
uated with the high-energy physics frequently used 
event generator pythia/jetset [195|. 



Whole spatial and spectral source distribution 

Once the spectrum of emitted antiprotons is known 
for a single PBH of given mass, the source term used 
for propagation is given by 



dEdtdV 



d^N, 



dEdt 



d^n 
dMdV 



dM , (52) 



where d^n/dMdV is the PBH mass spectrum to- 
day. It can be deduced from the initial mass spec- 



As shown in the main uncertainties on the esti- 
mation of the antiproton fiux from PBHs are associ- 
ated with astrophysical parameters. Several possible 
effects due to the quite large size of the horizon mass 
at the end of the inflation period (mostly imposed 
by gravitinos and moduli fields constraints) and to 
the possible photosphere that could form around hot 
PBHs were studied in and show that the present 
estimates could be substantially revised in the future. 
The emitted antiproton spectrum is mostly resulting 
from PBHs with masses between lO^^ g and 10^* g: 
the lightest ones are not numerous enough because 
of the shape of the nowadays mass spectrum and 
the heaviest ones do not emit much because of their 
low temperature. 



Density exclusion criterion Superimposed with 
BESS, CAPRICE and AMS data, the full antiproton 
flux, including the secondary component and the pri- 
mary component, is shown in Fig. ^ for 20 values 
of P0^^ logarithmically spaced between 5 x 10~35 
and 10^32 gcm^3 with fixed astrophysical parame- 
ters. The lowest curves are clearly in agreement with 
data whereas the upper ones contradict experimen- 
tal results. To derive a reliable upper limit, and to 
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antiproton kinetic energy (GeV) diffusion halo thickness L (kpc) 



Figure 33: Experimental data from BESS 95 (filled 
circles), BESS 98 (circles), caprice (triangles) and 
AMS (squares) superimposed with mean theoretical 
PBH spectra for between 5 • 10~^^ g cm~^ 

(lower curve) and 10~^^ g cm~^ (upper curves). 



a generalized as 



X 



^ (ar+^*+(Q«)) 

($'h(g^) _ $0)2 



(ar+cT"-(Qoy 

where cr*''+ and (7°+ ((t*''~ and cr'^^) arc the theoret- 
ical and experimental positive (negative) uncertain- 
ties. An upper limit on the primary flux for each 
value of the magnetic halo thickness is computed. 
The theoretical errors included in the function 
come from nuclear physics {p + He p + X and 
He + He p + X jSSj), and from the astrophysi- 
cal parameters analysis which were added linearly, in 
order to remain conservative. The resulting x^ leads 
to very safe results as it assumes that limits on the 
parameters correspond to 1 sigma. 



Resulting constraints Fig. ^ gives the upper 
limits on the local density of PBHs as a function of 
L. It is a decreasing function of the halo size because 
a larger diffusion region means a higher number of 
PBHs inside the magnetic zone for a given local den- 
sity. Between L = \ and L = 15 kpc (extreme astro- 
physical values), the 99% confidence level upper limit 
goes from 5.3 • 10"^"^ g cm~^ to 5.1 • 10""^^ g cm""^. 

Such limits are of great cosmological relevance. 
They are a unique probe of the early Universe on 
very small scales that could be equaled only, in a 
far future, by the detection of primordial gravita- 



Figure 34: Upper limits on the local density of PBHs 
as a function of the magnetic halo thickness L. 

equivalent, limits coming from the 100 MeV gamma- 
ray background, numerous constraints were derived 
on t he s pectral index of scalar perturbations (see, 
e.g. 1193): n < 1.27. A too "blue" spectrum would 
lead to an overproduction of PBHs in conflict with 
observations. 

This upper limit on the local PBH density ob- 
tained with antiprotons can also be used to severely 
constrain models with Broken Scale Invariance 
(BSI). A jump in the first derivative of the inflaton 
potential [198| should lead to a huge increase of the 
PBHs formation probability when the correspond- 
ing scale re-enters the horizon [199 . The antiproton 



limit can, therefore, be directly translated into an 
exclusion area in the parameter space of BSI infla- 
tionary models |200|. Finally, it allows interesting 



prospects for PBH dark matter [201| or Planck relics 
investigations |p02|. 



7.4.3 A new window for detection: antideuterons 

To go beyond an upper limit and try to detect PBHs 
it seems very interesting to look for antideuterons. 
Below a few GeV, there is nearly no background for 
kinematical reasons [162| and the possible signal due 
to PBHs evaporation could be easy to detect. The 
emission scheme is nearly the same but the proba- 
bility that an antiproton and an antincutron merge 
into an antideutcron is taken into account. The pos- 
sible detection range for the AMS experiment [ p03[ 
can be evaluated. It is shown on Fig. ^ func- 
tion of the three unknown parameters: L, the height 
of the magnetic halo, the coalescence momentum 
to form an antideuteron, and Pq^^, the local density 
of primordial black holes. The sensitivity of the ex- 
periment should allow, for averaged parameters, an 
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tor of six, if not a positive detection. The situation 
is very different than for antiproton: the main Umi- 
tation is due to the experimental sensitivity and not 
to the unavoidable physical background. Great im- 
provements can, therefore, be expected in the future 
and this investigation means seems very promising. 
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Figure 35: Parameter space (halo thickness L: 1- 
15 kpc; coalescence momentum pn- 60-285 MeV/c; 
PBH density p©: 10^^^ — 10^^^ g cm~^) within the 
AMS sensitivity (3 years of data taking) . The allowed 
region lies below the surface. 



8 Summary, conclusions and perspec- 
tives 

A consistent framework to understand the propa- 
gation of CR nuclei in the energy range 100 MeV~ 
100 GeV was presented in this paper. The observed 
fluxes of most species can be explained by assuming 
that once emitted from some sources located in the 
galactic disc, these nuclei undergo a diffusive prop- 
agation altered by escape through the boundaries, 
spallations, reacceleration, energy losses and galac- 
tic wind. The magnitude of these effects has been 
constrained using the B/C data, and the consistency 
of the model has been tested against the observed an- 
tiproton flux and by the study of radioactive species. 

This well-tested model has then been used to study 
the propagation of cosmic rays of a more hypotheti- 
cal origin, such as light antinuclei produced by SUSY 
galactic Dark Matter or Primordial Black Holes. In 
particular, limits on the abundance of primordial 
black holes in the galactic halo could be set. 

There are several directions in which this work 
may be extended. First, the constraints on the prop- 
agation parameters could be refined by considering 
other species, stable or secondary. However, this ap- 
proach is currently limited by the accuracy of the 
available data on cosmic ray fluxes and on the nu- 
clear cross sections. Second, a specific study of the 
EC unstable species could provide valuable informa- 
tion about the processes responsible for the acceler- 



the SUSY induced antiproton signal can be bettered 
by using the constraints on the propagation param- 
eters in a fully consistent way. Finally, the propaga- 
tion code we use should ultimately be able to yield 
the flux of all cosmic ray species (including gamma 
rays, electrons and positrons) at every position in the 
Galaxy. 



Local properties 
in our Galaxy 



Propagation 
in the Galaxy 



(8) Spatial Origin 



Astro-particle component 

^„ « p prim. SUSY '-. 

J" """ "" "-- ---«■ [ p prim.PBH\ 4) 
--- prim^PBH] (7f 
'^f d ;«im. SUSY :*** 

(-1) 

c: Extragalactic ? : 



Extension to PeV energy, 
CRs in the whole Galaxy 



: Diffuse 7 rays --- 
: Local 



■ The knee? 
^ {1 PEV) 



Figure 36: Schematic view of the subjects discussed 
in this paper. The stars indicate collaborations of 
LAPTH (Annecy-le-Vieux, France) members with ISN 
(★, Grenoble, France), lAP (★★, Paris, France) or infn 
(★ ★ Turin, Italy). Dashed boxes represent future 
projects. The numbers in parenthesis represent the 
publications. The starting point (1) is the first use 
of the elaborate propagation m odel presented here 
l4Sj, ( 2) is @ (3) is II, (4) is (lllL (5) is §^,(6) 
is U, (7) is H and (8) is l65|, |186| . The pre-|| 
works a re la belled as (-1) for |152| (-2) for and 
(-3) for jlTl. 
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